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Abstract 

Two-nucleon axial charge and current operators are derived in chiral effective field theory up to 
one loop. The derivation is based on time-ordered perturbation theory, and accounts for cancel¬ 
lations between the contributions of irreducible diagrams and the contributions due to non-static 
corrections from energy denominators of reducible diagrams. Ultraviolet divergencies associated 
with the loop corrections are isolated in dimensional regularization. The resulting axial current is 
finite and conserved in the chiral limit, while the axial charge requires renormalization. A com¬ 
plete set of contact terms for the axial charge up to the relevant order in the power counting is 
constructed. 

PACS numbers: 21.45.-v, 23.40-s 
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I. INTRODUCTION 


Chiral symmetry is an approximate symmetry of Quantum Chromodynamics (QCD), 
the fundamental theory that describes the interactions of quarks and gluons—the symmetry 
becomes exact in the limit of vanishing quark masses. Chiral effective held theory (yEFT) is 
the theoretical framework that permits the derivation of nuclear potentials and electroweak 
currents from the symmetries of QCD—the exact Lorentz, parity, and time-reversal symme¬ 
tries, and the approximate chiral symmetry. Pions and nucleons (and low-energy excitations 
of the nucleon, such as the A isobar), rather than quarks and gluons, are the degrees of free¬ 
dom of yEFT. Chiral symmetry requires the pion to couple to these baryons, as well as 
to other pions, by powers of its momentum Q and, as a consequence, the Lagrangian de¬ 
scribing their interactions can be expanded in powers of Q/A^, where ~ 1 GeV is the 
chiral symmetry breaking scale. Classes of Lagrangians emerge, each characterized by a 
given power of Q/A^, or equivalently a given order in the derivatives of the pion held and/or 
pion mass factors, and each containing a certain number of unknown parameters, the so 
called low-energy constants (LECs). These LECs could in principle be calculated from the 
underlying QCD theory of quarks and gluons, but the non-perturbative nature of this theory 
at low energies makes this task extremely difficult. Hence, in practice, the LECs are hxed 
by comparison with experimental data, and therefore effectively encode short-range physics 
and the effects of baryon resonances, such as the A isobar, and heavy-meson exchanges, not 
explicitly retained in the chiral Lagrangians. 

Within yEFT a variety of studies have been carried out in the strong-interaction sector 
dealing with the derivation of two- and three-nucleon potentials [1-9] and accompanying 
isospin-symmetry-breaking corrections [10-13], and in the electroweak sector dealing with 
the derivation of parity-violating two-nucleon potentials induced by hadronic weak interac¬ 
tions [14-17] and the construction of nuclear electroweak currents [18-25]. Most of these 
studies have been based on a formulation of yEFT in which nucleons and pions are the 
explicit degrees of freedom. A few, however, have also retained A isobars as explicit degrees 
of freedom. 

In this paper, the focus is on nuclear axial charge and current operators. These were 
originally derived up to one loop in heavy-baryon covariant perturbation theory (HBPT) in 
a pioneering work by Park et al. [18]. Here we re-derive them by employing a formulation 
of time-ordered perturbation theory (TOPT), which accounts for cancellations occurring 
at a given order in the power counting between the contributions of irreducible diagrams 
and the contributions due to non-static corrections from energy denominators of reducible 
diagrams [20]. Because of the different treatment of reducible diagrams in the HBPT and 
TOPT approaches, we hnd differences between the operators obtained in these two for¬ 
malisms as well as additional differences due to the omission of a number of contributions 
in Ref. [18], as discussed in Sec. VH. 

An accurate theory of nuclear electroweak structure and dynamics is relevant in several 
areas of current interest. One such area is that of low-energy tests of physics beyond the 
Standard Model in /d-decay experiments [26]. Phenomenologically, the weak interactions are 
known to couple only to left-handed neutrinos, and to violate parity maximally. However, 
beyond the Standard Model (BSM) theories have been constructed in which small deviations 
from these properties are introduced. These deviations affect the correlation coefficients 
entering /d-decay rates, and can in principle be detected. For a proper interpretation of 
these measurements and, in particular, to unravel possible signatures of BSM physics, it is 
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crucial to have control of the nuclear structure and weak interactions in nuclei. 

Another area of interest is that of neutrino interactions with nuclei and neutron matter. 
The low-energy inelastic neutrino scattering from nuclei is important in astrophysics and for 
neutrino detectors. The spallation of neutrons from nuclei by neutrino interactions is relevant 
in setting the neutron to seed ratio in core-collapse supernovae. Accurate predictions for 
neutrino-nucleus scattering cross sections, specihcally from the argon nucleus, are key to the 
measurements of supernovae neutrino fluxes, a major component of the Deep Underground 
Neutrino Experiment (DUNE). At temperatures of a few MeV, neutrino processes are also 
very important in core-collapse supernovae. One signihcant issue is the decoupling of various 
flavors of neutrinos and antineutrinos at the surface of the proto-neutron star. This sets 
the initial temperatures (flux versus energy) of e, /i and r neutrinos and antineutrinos. 
Understanding this initial flux is critical to interpreting the subsequent evolution of neutrinos 
and their role in the r-process. Neutrino and antineutrino interactions in neutron matter are 
also of importance in understanding the evolution of the very neutron-rich matter formed 
in neutron-star mergers, since they can potentially alter the neutron to proton ratio and 
signihcantly impact the r-process in neutron star mergers, currently considered to be an 
important source for r-process nucleosynthesis. 

The present paper is organized as follows. In Sec. II pion-nucleon [ttN) and pion-pion 
(tttt) interaction Hamiltonians are constructed from the chiral Lagrangian formulation of 
Refs. [27, 28]—for convenience these Lagrangians are listed in Appendix A, where a number 
of details relative to the construction of the Hamiltonians up to the relevant chiral order 
are also provided. In Sec. HI the power counting scheme and TOPT formulation adopted 
in the present work are described. These along with the interaction vertices obtained in 
Appendix B are utilized to derive two-nucleon axial charge and current operators up to 
one loop in Secs. IV and V, respectively. Ultraviolet divergencies associated with the loop 
corrections are isolated in dimensional regularization: the resulting axial current is then 
found to be hnite, while the axial charge requires renormalization. All this along with the 
renormalization of the one-pion-exchange (tree-level) axial charge is discussed- in Sec. VI. A 
number of details are relegated to Appendix C, where a complete set of contact terms for the 
axial charge (up to the relevant order) is constructed, to Appendix D, where loop functions 
entering the axial current are dehned, and to Appendix E, where a listing of counter-terms 
is given. In Sec. VH a summary and discussion of our results as well as a comparison 
between the expressions for the axial operators obtained here and those of Park et al. [18] 
are provided. Conclusions are summarized in Sec. VHI. 


II. INTERACTION HAMILTONIANS FROM CHIRAL LAGRANGIANS 


The chiral Lagrangian describing the interactions of pions and nucleons is given by 


C 


C 


ttN 


r 


( 2 . 1 ) 


where 

HttW ^ + cfh + cfl + ... , (2.2) 

£„ = £S + 4t + ... . (2.3) 

and the superscript n specihes the chiral order {Q denotes generically the low-momentum 
scale), i.e., the number of derivatives of the pion held and/or insertions of the pion mass. 
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External fields are counted as being of order Q. Since we are interested in deriving nuclear 
potentials and currents up one loop, it suffices to retain in C up to and CnJ. The 
Lagrangians (in fact up to order n = 4) and have been given, for example, in 
Refs. [27] and [28], respectively, and are listed in Appendix A of the present paper for 
completeness. The total Lagrangian can be written as 

C = N {i$ - m + T^^ doTTa + K + A) A^ 

+ ^ Gab OoTTb + <9Va Gab OiTlb “ ml TlaHab TTb^ “ U Fab {O^^TTb) , (2.4) 

where tTq is the pion held of isospin component a, N is the iso-doublet of nucleon helds, 
is the axial-vector held of isospin component a, is the pion decay constant, and m 
and ttItt are, respectively, the nucleon and pion masses. The symbols r°, A[j, and A denote 
combinations of the pion and axial-vector helds (and their derivatives) and/or of pion mass 
factors, having the following expansions 

r; = r;:(o) + r;(i) + r;(2), (2.5) 


and similarly for A)^, and 

A = A(l)+A(2) + A(3) , (2.6) 

where the argument n in r[](n), A[,(n), and A(n) specihes the power counting Q”. The 
symbols Gab, Gab, Hab, and Fab denote three-by-three matrices in isospin space, containing 
powers of the pion held and/or pion mass. A listing of all these quantities, limited to the 
terms relevant for the construction of the currents at one loop, is provided in Appendix A. 
At this stage the various helds, masses, and coupling constants are to be understood as bare 
(un-renormalized) quantities. 

From the Lagrangian C in Eq. (2.4) the conjugate momenta relative to the pion and 
nucleon helds follow as 




dC 

d{doN) 

dC 

didoTla) 


Gabd%b-UFabAl + NTlN , 


and the Hamiltonian then reads 


( 2 . 7 ) 

( 2 . 8 ) 


w = n* d„N + n„ SoTT. 


( 2 . 9 ) 


where "Ho, 

■Ho = ^ (Ha Ua - d^TTa dill a + m^ T^a T^a) + N {-% Oi + m) N , (2.10) 

is the free pion and nucleon Hamiltonian, while is the Hamiltonian accounting for the 
interactions between pions and nucleons as well as between these and the external held. By 
only keeping terms linear in the latter, the interaction Hamiltonian is given by 

w, = )n„ [(G-')„, - s^] n, - i [n. (Vr» n) + h.c.] 

+y [n. (G-‘)„,F,,A? + h.c.] - h [(ivrSiv) (G-')„,F^2i; + h.c.] 
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( 2 . 11 ) 


+i (/V r; N) (G-)„J {N r'>N)-N (A- 8,7,, + A) N 

\ \ 7Tl^ 

It admits the following expansion in powers of Q\ 

Uj = nf^+Hf+Hf + ... , ( 2 . 12 ) 

and the vertices corresponding to the various interaction terms are listed in Appendix B. 


III. FROM AMPLITUDES TO CURRENTS 

The expansion of the transition amplitude for a given process is based on TOPT. Terms 
in this expansion are conveniently represented by diagrams. We distinguish between re¬ 
ducible diagrams (diagrams which involve at least one pure nucleonic intermediate state) 
and irreducible diagrams (diagrams which include pionic and nucleonic intermediate states). 
The former are enhanced with respect to the latter by a factor of Q for each pure nucleonic 
intermediate state (see below). In the static limit—in the limit m —)■ cx), i.e., neglecting 
nucleon kinetic energies—reducible contributions are infrared-divergent. The prescription 
proposed by Weinberg [29] to treat these is to dehne the nuclear potential and currents as 
given by the irreducible contributions only. Reducible contributions, instead, are generated 
by solving the Lippmann-Schwinger (or Schrodinger) equation iteratively with the nuclear 
potential (and currents) arising from irreducible amplitudes. 

The formalism developed by some of the present authors is based on this prescription [20]. 
However, the omission of reducible contributions from the dehnition of nuclear operators 
needs to be dealt with care when the irreducible amplitude is evaluated under an approxi¬ 
mation. It is usually the case that the irreducible amplitude is evaluated in the static limit 
approximation. The iterative process will then generate only that part of the reducible 
amplitude including the approximate static nuclear operators. The reducible part obtained 
beyond the static limit approximation needs to be incorporated order by order—along with 
the irreducible amplitude—in the dehnition of nuclear operators. This scheme in combina¬ 
tion with TOPT, which is best suited to separate the reducible content from the irreducible 
one, has been implemented in Refs. [21, 23, 25] and is briehy described below. The method 
leads to nuclear operators which are not uniquely dehned due to the non-uniqueness of the 
transition amplitude off-the-energy shell. While non unique, the resulting operators are 
nevertheless unitarily equivalent, and therefore the description of physical systems is not 
affected by this ambiguity [23, 30]. 

We note that an alternative approach, implemented to face the difficulties posed by the 
reducible amplitudes, has been introduced by Epelbaum and collaborators [31]. The method, 
referred to as the unitary transformation method, is based on TOPT and exploits the Okubo 
(unitary) transformation [32] to decouple the Fock space of pions and nucleons into two sub¬ 
spaces, one containing only pure nucleonic states and the other involving states that retain 
at least one pion. In this decoupled space, the amplitude does not involve enhanced con¬ 
tributions associated with the reducible diagrams. The subspaces are not-uniquely dehned, 
since it is always possible to perform additional unitary transformations onto them, with 
a consequent change in the formal dehnition of the resulting nuclear operators. This, of 
course, does not ahect physical representations. 
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The two TOPT-based methods outlined above lead to formally equivalent operator struc¬ 
tures for the nuclear potential and electromagnetic charge and current up to loop corrections 
included, which makes it plausible to conjecture that the two methods are closely related. 
However, this topic has not been investigated further. In what follows, we focus on the 
method developed in Refs. [21, 23, 25] and show how nuclear operators are derived from 
transition amplitudes. Here, we are especially interested in the construction of the two-body 
weak axial charge and current operators. We will not discuss the aforementioned unitary 
equivalence between operators corresponding to different off-the-energy-shell extrapolations 
of the transition amplitudes. This issue has already been addressed in considerable detail 
in Ref. [23] for the case of the two-body nuclear potential and electromagnetic charge and 
current operators. Similar considerations apply to the present case. 

The starting point is the conventional perturbative expansion for the amplitude 

oo / .. \ n—1 

where | i) and | /) represent the initial and hnal states, respectively | N 1 N 2 A) and | N[N!^ 
{A denotes generically the external axial held), of energies Ei and Ef with Ei = Ef, Hq is 
the Hamiltonian describing free pions and nucleons, and Hj is the Hamiltonian describing 
interactions among these particles {Hq = f dx'Ho(x) and similarly for Hj, with Hq and H/ as 
dehned in Sec. H with the various helds taken in the Schrodinger picture). The evaluation of 
this amplitude is carried out in practice by inserting complete sets of Hq eigenstates between 
successive terms of H/. Power counting is then used to organize the expansion in powers of 

Q/A^ < 1 . 

In the perturbative series, Eq. (3.1), a generic (reducible or irreducible) contribution 
is characterized by a certain number, say N, of vertices, each scaling as x 
(f=l,..., iV), where ctj is the power counting implied by the specihc term in the inter¬ 
action Hamiltonian Hj under consideration and /dj is the number of pions in and/or out of 
the vertex, a corresponding iV-1 number of energy denominators, and L loops. Out of these 
A^-1 energy denominators, Nk of them will involve only nucleon kinetic energies and possi¬ 
bly, depending on the particular time ordering under consideration, the energy Uq associated 
with the external held, both of which scale as while the remaining N — — 1 energy 

denominators will involve, in addition, pion energies, which are of order Q. Loops, on the 
other hand, contribute a factor each, since they imply integrations over intermediate 
three momenta. Hence the power counting associated with such a contribution is 

N 

J-J gai-ft/2 
i=l 

Clearly, each of the N — Nk — 1 energy denominators can be further expanded as 

1 ^ [1 + (^^ - Elf ^ 

Ei Ej UJ'jY UJ-jz _ bJ'jx 

where Ej denotes the energy of the intermediate state (including the kinetic energies of the 
two nucleons and, where appropriate, the energy of the external held), and the pion 
energy (or energies, as the case may be)—the ratio {Ei — Ei)/un is of order Q. The leading 
order term —I/utt represents the static limit, while the sub-leading terms involving powers 



X [qAn-Nk-i) g-2AfKj X g3L 
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of {Ei — Ej)/ut^ represent non-static corrections of increasing order; elsewhere [ 20 , 21 ], we 
have referred to these as recoil corrections. 

Interactions with the external axial held are treated in hrst order in Eq. (3.1), and in¬ 
spection of the Q scaling of the various terms shows that the associated transition amplitude 
admits the following expansion 

= + ... , (3.4) 

where is of order Q^. Next, we denote the two-nucleon strong-interaction potential 
with V and the weak-interaction potential with ^5 = ■ js.a, where and 

are, respectively, the nuclear weak axial charge and current operators and = (A°, A^) 
is the external axial held. We construct n -|- ^5 by requiring that iterations of n -|- us in the 
Lippmann-Schwinger equation [23] 

(n -f Us) -f (n + Us) Go {v + '^ 5 ) + ("^ + v^) Gq {v + ^ 5 ) Gq {v + v^) + ... , (3.5) 


match the T 5 amplitude, on the energy shell Ei = Ef, order by order in the power counting; 
here Go denotes the propagator Go = l/(E'i — Ej + ip). The potentials v and V 5 have the 
following expansions 


Us = -F . . . 


(3.6) 

(3.7) 

where the potentials have been derived in Refs. [21, 23], in particular vanishes [23], 


and ■ jg”j. The superscript (n) on V 5 and Tg only refers to the power 

counting of pg*)] and jg^j, and does not include the power of Q associated with the external 


held. The matching between Tg""^ and Ug"''* leads to the following relations [23] 


(n) 


(-3) 

-'s 

,(- 2 ) 

25 


-i(—3) 
5 

-i(-2) 


Go Go 4“^^ 


G — -^5 


,,(o) _ 7 ^( 0 ) _ 

^5 — -^5 


4 Go Go -I- permutations 

4"^^ Go Go 4“^^ , 

4~^^ Go Go Go -I- permutations 
4~^^ Go Go -I- permutations 

4“^^ Gon^°^ Go 4”^^ 

4”^^ Go Go 4”^^ 


(3.8) 

(3.9) 

(3.10) 


(3.11) 



7 ^( 1 ) 

-'s 


4 ~^^ Go Go Go Gq + permutations 

4 ~^^ Go Go Go + permutations 

4 Go Go + permutations 
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— Go Go + permutations 

- 4”^^ Go Go 4”^^ 

- 4"^^ Go Go 4"^^ , (3-12) 

and a similar set of relations is obtained between TG) and v^"'\ i.e., the amplitudes and po¬ 
tentials in the presence of strong interactions only [23]. These relations allow us to construct 
t;G) and 4”^ from TG) and 



a1 a2 c 


FIG. 1. Diagrams al and a2 contribute to the one-body axial current operator at order QGs) 

Diagram c contributes to the one-body axial charge operator at order QG"^) , Nucleons, pions, 

and axial fields are denoted by solid, dashed, and wavy lines, respectively. Only a single time 

( 2 ) 

ordering is shown for diagrams a2 and c. The full dot in c is from the interaction vertex 
see Appendix B. 

The weak axial charge and current operators at leading order consist of the single-nucleon 
contributions shown in Fig. 1 and are given by 

Pv Hq) = D,acri ■ (p'l + pi) (27r)^(5(pi + q - p[) + (1^2) , (3.13) 

’ 4m 

Jv ^(q) = -^ D,a <Ti - ^ 2 • q (27r)^5(pi + q-p'i) + (1^2) , (3.14) 

’ 2 I 

where q is the momentum carried by the external held, and pj and p) are the initial and 
hnal momenta of nucleon i. The counting Q~^ of (panel al in Fig. 1) follows from the 
product of a factor associated with the ANN current vertex (recall that the Q scaling 
of the external held is not counted), and a factor Q~^ due to the momentum-conserving <5- 
function 5 (p 2 — P 2 ) implicit in disconnected terms of this type. Evaluation of the pion-pole 
contribution (panel c), in which the axial source couples directly to the pion which is then 
absorbed by the nucleon, leads to the expression in Eq. (3.13). In this disconnected 
term, the counting Q~'^ accounts for the Q~^ factor due to 5 (p 2 — P2), the factors Q and 
of the ttA and ttNN vertices, respectively, and the factor Q~^ from the pion held normal¬ 
ization and energy denominator associated with the intermediate state. A similar counting 
is applied to panel a2 in Fig. 1 contributing to js^a. 

There is no direct coupling of the nucleon to the interaction —(p^/2)A^ r ■ Aq 7 ^ 7 ^ N 
in 

-NA{2)N , 
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with A( 2 ) as given by Eq. (A59) occurs with the opposite sign in 

-(/,/ 2 ) [jvr;(l)/v + , 

with r°(l) as in the hrst term of Eq. (A57) and {G)~l = Fat = Sab up to Tia'^b or terms, 
and hence cancels out in Eq. (2.11). The single-nucleon axial charge of the correct sign and 
strength follows from the sum of the two time-ordered contributions of diagram c with the 
full dot representing the interaction ( 5 f^/ 2 /^) N r ■ 11 7 ^ 7 ^ N from 

-(1/2) [n4G-')^j Afr»(i)Af + h.c.] . 


Because of the different power counting of the leading order terms in the current and 
charge operators, the strong interaction potentials needed to construct these operators up 
to order n = 1 include corrections up to n = 3, i.e., in the case of the current and up 
to n = 2, i.e., in the case of the charge. The leading order (LO) term consists of 
(static) one-pion-exchange (OPE) and contact interactions, while the next-to-leading order 
(NLO) term (as already noted) vanishes (see Ref. [23]). The next-to-next-to leading 
order (N2LO) term contains two-pion-exchange (TPE) and contact interactions, the 
latter involving two gradients of the nucleon helds. The term was originally derived in 
Ref. [1], and is well known. However, at N2LO there is also a recoil correction to the OPE, 
which we write as [30] 


7 ^)(„) = „(»)(k) EME '2 - E,) 


27 


(3.15) 


where nl°^(k) is the leading order OPE potential, dehned as 


n(°)(k) = 


9 a h h ^ 

Ti ■ T2 CTi ■ k 0-2 ■ k ^ , 


4 /. 


U)t 


(3.16) 


Ei (pi) and E' (p') are the initial and hnal energies (momenta) of nucleon i, and k = pi — Pi. 
There is an inhnite class of corrections vir (u), labeled by the parameter z/, which, while 
equivalent on the energy shell {E[ + E 2 = Ei + E 2 ) and hence independent of u, are different 
off the energy shell. Friar [30] has in fact shown that these different off-the-energy-shell 
extrapolations are unitarily equivalent. 

The next-to-next-to-next-to-leading order (N3LO) term includes interactions gen- 
erated by vertices from the sub-leading Lagrangian —these are of no interest for the 

present discussion—as well as non-static corrections to the N2LO potentials v^'^\ Among 

/o\ 

these, the TPE correction ^271 (^) (froui direct and crossed box diagrams) depends on the 

/Q\ 

specihc choice made for Vn (i^)- However, as shown in Ref. [23], the unitary equivalence 
remains valid also for In the derivation of the axial current at n = 1 below, the 

choice z/ = 0 is made for and specihcally Eq. (3.15) above and Eq. (19) of 

Ref. [23]. The remaining non-static corrections in the potential are as given in Eqs. (B 8 ), 
(BIO), and (B12) of that work. Clearly, different choices in the off-the-energy-shell extrap¬ 
olations of these potentials will lead to different forms for (some of) the j 5 ^a(z^) corrections 
to the axial current. As shown in the case of the electromagnetic charge operator [23], one 
would expect these different forms to be unitarily equivalent. However, this has not been 
verihed explicitly in the present case. 
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IV. AXIAL CHARGE 


The nuclear weak axial charge two-body operator can be written as 


^ _ OPE , TPE , OT 

Pb,a — P5^a + P5,a + P5,a 


(4.1) 


namely as a sum of terms due to one-pion exchange (OPE), two-pion exchange (TPE), and 
contact contributions (CT). We defer the discussion of loop corrections to the OPE axial 
charge (and current) and of their renormalization to a later section. In the following, and in 
Sec. V as well, contributions to the OPE and TPE (or MPE in Sec. V) operators are labeled 
by the power counting superscript (n). While each individual contribution is not explicitly 
identihed as being OPE and TPE (or MPE), this is obvious from the context. 

Here and throughout this paper, we adopt the following conventions. The momenta k* 
and Kj are dehned as 

Ki = (p' + Pi) /2 , ki = p' - Pi , (4.2) 

where pi (p') is the initial (hnal) momentum of nucleon i. A symmetrization (1 ^ 2) and 
an overall momentum-conserving (5-function (27r)^(5(ki -|- k 2 — q) are understood in all terms 
listed below unless otherwise noted. 


A. Leading one-pion and two-pion exchange contributions 


OPP TPP 

Diagrams contributing to pYf leading order and to p^f are shown, respectively, 
in panels al and a2, and panels cl-cl2 of Fig. 2. The contributions of al-a2, and cl-c2 and 
c5-c6 are given by 


Pv Vl) 

pS(c1 + c2 ) 
pS(c 5 + c6) 


■Pa . . ,1 

^ T 2 ),cr 2 -k 2 ^ , 

Pv^al) , 

* ^ 'r2)aCri ■ k2l^°\k2) , 

* 4ri,aCTii (0-2 X k2)^. jff\k2) 


+ (Ti X T2)a [kj J^^\k2) - ^^^^^ 2 )] CTi ■ k2 


(4.3) 

(4.4) 

(4.5) 


(4.6) 


while those of c3-c4, c7-c8, and c9-cl2 vanish. Corrections proportional to 1/m to topologies 
al and a2, due to non-static corrections to the energy denominators, that enter at order Q, 
vanish after summing over all time orderings. Contributions, coming from and P-L^^tInn^ 

to topologies al and a2, that enter at order Q, turn out to vanish. The freedom in the choice 
of pion held, parametrized by the parameter a in Appendix A, introduces an (^-dependence 
in the interaction vertices with three or four pions, see Appendix B. The contributions of 
diagrams c4 and c8, which include a 37r vertex, turn out to vanish identically. But in general 
this a dependence must cancel out exactly, as is indeed the case for the two-nucleon axial 
charge and current operators obtained in this work. The loop functions have been dehned 
as 

= (4-7) 
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J(°)(A;) = 

/ (222)3 

(4.8) 

J^^\k) = 

/ (2,2)3 

(4.9) 

= 


(4.10) 


with 


f{u_,u+) = 
g{uj-,uj+) = 


0 ;+ a;_ (a;+ + a;.) ’ 
OJ-y. OJ— -\- Upi_ 


and 


0;+(x;5?.(a;++ a;_) 
w± ^ ^(p ±k )2 +4ml . 


(4.11) 

(4.12) 

(4.13) 


They are evaluated in dimensional regularization [21], Insertion of the hnite parts of these 
loop functions leads to 


Pv(cl + c 2 ) = -i 
P^lic5 + c 6 ) = -i 


9A 


128 7r2/4 


. X , S2 , fs2 + k2 

Ti X T 2 )aCri ■ k 2 — In 


9a 


128 ttV." 

Tl X T2)aCri ■ k2 


k2 \S2 - h 

A ! X 1 -52 , <§2 + /i (2 

4ri,a {(Tl X <72) ■ k2 — In 


(4.14) 


^2 S2 ~ ^2 
/c| + 2 §2 j $2 ~\~ k2 
^2 S2 S2 ~ ^2 


where 


The divergent parts read 


p£(c 1 + c 2 )|oo = -i 


i = A/4m2+k 2 . 


TT ' ] 


9A 


p£(c5 + c 6 )|oo = -i 


1287r2/4 

9a 


[ti X T2)aCri ■ k2 - 1) , 


(4.15) 


(4.16) 


(4.17) 


327r2/, 


4 

TT L 


Tl,a {(Tl X <72) ■ k 2 d 


3 /I 

I (’’■ ^ .^ 2)«'^1 ■ k 2 I 4< + g 


(4.18) 


with the constant d. dehned as 


9 771^ 

d. = —+7-ln47r + ln^-l , 
e 


(4.19) 


where e = 3 — d (d is the number of dimensions), 7 is Euler’s constant, and /i is a renormal¬ 
ization scale. 
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FIG. 2. Diagrams contributing to the OPE axial charge at leading order Q~^ (panels al and a2), 
and to the TPE axial charge operator at order Q. Nucleons, pions, and axial fields are denoted by 
solid, dashed, and wavy lines, respectively. Only a single time ordering is shown for each topology. 


B. Contact contributions 

At order there are no contact terms contributing to Those at order Q are given 

by (see Appendix C) 

4 

= (4.20) 

i=l 

where the Zi are (unknown) LECs and the operators Oj with i = 1,... ,4, symmetrized with 
respect to the exchange 1^2, have been dehned as 

Oi = i (ti X T2)„ (cTi ■ k2 - 0-2 ■ ki) , (4.21) 

O 2 = i (ti X T2)„ (cTi ■ ki - <72 ■ k2) , (4.22) 

O 3 = i (<Ti X < 72 ) ■ (rpa k 2 - 72, a ki) , (4.23) 

O 4 = (ri,a - 72,a) (^1 - ^ 2 ) ' (Ki + K 2 ) . (4.24) 
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We observe that the loop divergencies from cl-c2 and c5-c6 can be reabsorbed in the LECs 
zi and Z 3 . 

V. AXIAL CURRENT 

Before considering the two-body contributions, we note that at order Q~^ there are 
relativistic corrections to the one-body current represented in diagrams bl and b2 of Fig. 3, 
given by 

j|"^^(bl) = cTi-h ^ki X Ki - cTi ■ Ki Ki + jCTi ■ kiki , (5.1) 

’ [2 4 

jya^Hb2) = — ^ ^ .Ki ki-Kij , (5.2) 

’ + ml L ’ 2 J 

where b2 contains two contributions at order Q~^\ one is from the 1/m^ terms originating 
from the non-relativistic expansion of the ttNN interaction the other is due to the 

1/m terms in and the (leading) non-static corrections (proportional to 1/m) to energy 

denominators. The bl current has been found to give a signihcant contribution to the cross 
section for proton weak capture on ^He of interest in solar physics [33]. 



FIG. 3. Diagrams illustrating the relativistic corrections to the one-body axial current. Nucleons, 
pions, and axial fields are denoted by solid, dashed, and wavy lines, respectively. Only a single 
time ordering is shown for diagram b2. See text for further explanations. 


As for the charge, the two-body current is written as a sum of one-pion exchange (OPE), 
multi-pion exchange (MPE), and contact (CT) terms (notation and conventions are as in 
Sec. IV), 


. _ -OPE , -MPE , -CT 

J5,a Js^a J5,a J5,a 


We discuss here. It is well known [33] that a single contact term occurs at order Q^, 
which we choose as 


= Xo[(Ti X Ts)^ CTi X 0-2 


+ ml 


(ti X Ts)^ q ■ (cTi X 0 - 2 ) ] , 


(where the second term of Eq.( 5.4) is the pion-pole contribution) and none at order Q (see 
Appendix C). This term is due to the interaction NN and, as hrst pointed 

by the authors of Ref. [34], the LEG Zq is related to the LEG cd (in standard notation) 
entering the three-nucleon potential at leading order. The two LEGs cd and ce which 
fully characterize this potential have been recently constrained by reproducing the empirical 
value of the Gamow-Teller matrix element in tritium /3 decay and the binding energies of 
the trinucleons [35, 36]. 
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A. Leading one-pion and multi-pion exchange and short-range contributions 


Leading contributions to and are shown, respectively, in panels dl-d2, and 

panels el-e23 of Fig. 4. There are no contributions at order Q~^ from diagrams dl and 
d2: in dl the interaction contains no coupling to the held Aq, while in d2 the sum 

over the 6 time orderings, when leading order vertices from and are 

considered, vanishes. The hrst non-vanishing contributions enter at order Q^, and read 


= ^ '^ 2 ), 


9a 


■ Ki Ce + 1 
i - 

2m 4m 


1 

CTi X q + ( C 4 + -— 1 CTi X k 
4m 


cr2 ■ k2 ^ 


, 9A , ,1 

+ C3 r2,a k2 <72 ■ k2 —^ , 


(5.5) 




9A q 


[ 7 - 2 ,a (4 Cl m^ + 2 C 3 q ■ k 2 ) - C 4 (ti x T 2 )„ <Ti ■ (q x k 2 )] <T 2 ■ k 2 ^ 

OJn 


2 /2 g 2 + ^2 

2 ^ 2 ('^1 ^ '^2)a (2 Ki + i <Ti X ki) ■ (q + k 2 ) <72 ■ ks ^ 




■ 9A 


q 


8 m + m^ 


(n X T2)„ (Ki ■ ki + 2 K 2 ■ k2) <72 ■ k2 ^ . 


Un 


(5.6) 


For the diagrams contributing to only a single time ordering is displayed for 

each topology. It is understood that denominators involving pion energies in the reducible 
topologies of diagrams el-e2, e6-e7, e8-el0, el3-el4, e20-e21, and e22-e23 are expanded as 
in Eq. (3.3). The resulting contributions depend on the off-the-energy-shell prescription 
adopted for the non-static corrections to the OPE, TPE, and OPE-contact potentials [23]. 
Different prescriptions lead to different formal expressions for these corrections as well as 
the accompanying weak axial current operators, which, however, are expected to be re¬ 
lated to each other by unitary transformations. This unitary equivalence was discussed in 
considerable detail in Ref. [23], where it was explicitly verihed to hold in the case of the 
electromagnetic charge operator. Here we reiterate that the axial current operators derived 
below are obtained by adopting the z/ = 0 prescription for the non-static corrections to the 
OPE, TPE, and OPE-contact potentials, as given in Eq. (3.15) of the present work and in 
Eqs. (19), (B8), (BIO), and (B12) of Ref. [23]. We hnd that the contributions of diagrams 
e3, e6-e7, ell-el4, el8-el9, e22-e23 vanish, while those of the remaining diagrams are given 
by 



9a 


16/.^ 

q 


^ 2 , 


Rlf(k2) (Tij - k2 R^^\k2) <7i ■ k2 


,2'jib'll). 


9a 


16/) 

9a 




kl R^^\ki) - R^^\ki) 


^2 , 


q 


32 // g2 + ^2 


T 2 ,a 


kl R'^^\ki) - R^‘^\ki) [(10 a - 1) <72 ■ k2 + <72 ■ ki] 


(5.7) 

(5.8) 

(5.9) 


- (ti X T2)„ (kl) (<7i X kl). a2,j , 


(5.10) 
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FIG. 4. Diagrams contributing to the OPE axial current operator at order and to the MPE 
axial current at order Q. Nucleons, pions, and axial fields are denoted by solid, dashed, and wavy 
lines, respectively. Only a single time ordering is shown for each topology. 


jS(e8) = 


9a 


16/^ 




(cTi X k 2 ) X k 2 


+ 


kl S^^\k2) - S^^\k2 


CTi - 


kl S^^\k2) - S^^\k2 
klSff{k2)-Sff{k2) 


aij 


+ 2 D,a (cr2 X k2) X k2 S^‘^\k2) 


(5.11) 
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(5.12) 


jS(e9) = 

jK(elO) = 




q 2 

q 2 7^2 (2 T 2 ,a - ri,a) (/C 2 )] CTi ■ ks 

+ (ti X T2)„ (ki) (0-2 X ka). cii^- , 

r2,a(10aq-3ki+ k 2 ) klR^^\ki) - R^^\ki) 

Jn I 


jS(el5) = 


(5.13) 


-4 (n X T 2 )„ (ki) (cTi X ki 


0-2 ■ k2 


(5.14) 


jS(«16) = 574;'■2.1 727777? 2 (bml + 2 kl + k^ + q‘^) [*■?«*“>(*■,) - 

O^JtT y + ^TT . 

+ [A;^i?(°)(A;i) -i?(^)(A:i)] - 20 a {q^ + kj + 2ml) Ikj R^^'>(k^) - R^^\ki) 

^ (72 ■ k2 

+ 80aJi2 -^ 


jU(el7) = ^ 


jK(e20) = 
jS(e21) = 


^16% g2 + rnl (o-i X ki), (k2 + q)^. , 

^T-2,a ") "712^^^ , 

8/4 g2 + m2 a;2 

Ct Ti^a <714 <72 , 

^ J 77 


-^-^q-j^(e20) , 


g2 + m2 


(5.15) 


(5.16) 


(5.17) 

(5.18) 


where the constants J^n are as in Eq. (B2), and the loop functions R^^ have been dehned 


R^^\k) = 


-f{u+,u_) , 


(5.19) 




«'''(*) - / 7^ 


-PiPj/(a;+,a;_) , 


^P^/(x^+,x^-) , 


(5.20) 

(5.21) 


(5.22) 


7?l;\k) = 


^PiPjP^ , 


(5.23) 

(5.24) 


/(x;+,x;_) = 


cai ca/ 


(5.25) 
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(5.26) 


The loop functions are defined similarly with /(a;+,a;_) replaced by 


^(cn+,cn_) 


ooX cot 


1 d 
Adm^ 


f{uj+,uj-) . 


After dimensional regularization, we obtain 


R^^\k) 

R^^\k) 


dz 


IGyrio M{k,z) 


/ dz 
AtT In 


M{k,z)-- 


1 z-z 




k^ 


-I-n L "" 


6 ij M{k, z] 


12 M(/c, 2 : 

1 ■ 
4 M{k,z) ^ \ 


RW{k) = - [ dz 


TT 


M(k, zf -1 (. - M(k, z) e + n 




kj 


M{k,zf-^{z-zfM{k,z)k^ 


21 

L 


(z-zyM{k,z)-- 


1 z-z 




28 M{k,z] 


k^ 


kikj 


(5.27) 

(5.28) 

(5.29) 

(5.30) 

(5.31) 


where 

M{k,z) = \/zzkP' + ml , (5.32) 

and 

z = l-;2. (5.33) 

The regularized S\^\k) loop functions easily follow from Eq. (5.26). Inserting these relations 
into the equations above, and noting that the a dependence cancels out upon summing the 
contributions of diagrams e5, el5, el6, and el7, we obtain the expressions reported in 
Appendix D. No divergencies occur in these loop corrections at order Q, consistently with 
the fact that there are no contact terms in the axial current at this order. Contributions 
coming from proportional to dj’s, that enter through topologies dl and d2 turn out to 
vanish. 


VI. RENORMALIZATION OF THE ONE-PION EXCHANGE AXIAL CHARGE 

We now proceed to renormalize the order Q loop corrections to the OPE axial charge 
operator (as shown below, no renormalization at this order is needed for the loop corrections 
to the OPE axial current). We hrst construct the set of relevant counter-terms, and then 
carry out the renormalization of the nucleon and pion masses, held rescaling factors and 
Ztv, pion decay constant /^r, nucleon axial coupling constant qa, and, lastly, loop corrections 
to the OPE axial charge. We dehne 

7r„=V^<, iV = v^iV^ (6.1) 

where and denote, respectively, the renormalized pion and nucleon helds, and and 
Ztv are the corresponding held rescaling constants, assumed to have the following expansions 

Zt^ = 1 - 1 - 5Zt^ , dZ^^ ~ , ( 6 . 2 ) 
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Zm = 1 + 5Zm , 5Zm ~ ■ 

We also define the physical pion mass and nucleon mass rrf as 

= ^l + ~ , 

rrf = m + 6 m , 6 m ~ . 

As illustrated in Appendix E, the total Lagrangian can be written as 


C 


N'^{t0-m^ + rl' doK + Al' + A') 

+1 (a“< gIj s„< + d'K GU d.K - < Ki f 

_ _ 

+(5m + 6ZjvN" - m’’) A" + ^ 7r> 


r 

a 






(6.3) 

(6.4) 

(6.5) 


( 6 . 6 ) 


which is then expressed in terms of renormalized helds and masses, but bare coupling con¬ 
stants qa and f-j^ and LECs. This Lagrangian has essentially the same form as the bare 
one in Eq. (2.4) (the primed quantities are dehned in Appendix E), and leads to a similar 
interaction Hamiltonian as in Eq. (2.11), 


1 -Li = "H/^Eq. (2.11) with primed quantities and renormalized helds and masses] 

_ _ 

- 6 mN"N^ - 6 ZnN" {i^di - m^) A" - . (6.7) 

In addition to the vertices listed in Appendix B, this Hamiltonian generates vertices cor¬ 
responding to the set of counter-terms in Eqs. (E9)-(E15), explicit expressions for which 
follow from those in Appendix B. 


A. Field and mass renormalization 


The determination of the scaling factors = 1 + 6 Zt^ and Zj^j = l + 6 ZAr for the pion and 
nucleon helds, and the renormalization of the pion and nucleon masses have been discussed 
recently and in considerable detail in Ref. [17]. We only quote the results here: 


6 mi = 21 


m 


r 4 


m 


r2 


f2 
J TT 

6 m = — 


+ 


Cl - 




01 ; 


J 


12 


6Z^ = -2 


6Zm = 


m 


r2 


8 gi 

8 n 


li + 


lOa- 1 


Ji 


01 ; 


J 


13 


( 6 . 8 ) 

(6.9) 


where the constants Jmn are dehned in Eq. (B2). Only leading corrections are provided 
above, but for 6 m which also includes the sub-leading term of order proportional to J 12 . 
The sign for 6 m dihers from that in Ref. [17], since there m’’ = m — 6 m. 


B. Renormalization of the pion decay constant 

The relevant interaction Hamiltonians are 

Hif = uJd^{A^-d,7z^ + A^-W) , ( 6 . 10 ) 


18 



m J! = ^ / dx 


2 a) A* ■ tt" tt" ■ diTT^ - (2 a + 1) A* ■ tt" ■ tt*" 


^SttA ~ 7)~r / — ' n n ' Uin — ui t ' ui a 

^Jn J ^ 

+2 (a - 1/2) A° < nl + 2a A° « tt^-W + W- tt; 


Hif=ld^ ^^^A*-a,7r’--^/4-A*-a,7r^ + A°-n’-) , 


( 6 . 11 ) 

( 6 . 12 ) 


/Q\/ ( ^'\ / 

where ^s-Ja same as in Eqs. (B42) and (B46) but in terms of renormal¬ 
ized pion field and mass, while relative to Eq. (B43) includes counter-terms. The 

contributions illustrated in Fig. 5 read 

al = -ifn (k ■ Aa - a;A°) , (6.13) 


a2 = -;^ Joi [ - (5 a + 1/2) A, ■ k - (5 a - 3/2)A° J , 

a3 = -2i ^^ k ■ Aa -1- (-k ■ A^ - a;A°) . 

Jtt ^ 

We now require that the renormalized (physical) pion decay constant is equal to 


(6.13) 

(6.14) 


(6.15) 


a1 a2 


FIG. 5. Diagrams relevant for the renormalization of f-,^ 


implying 


— i/)) (k ■ A — 00 A^) = al -I- a2 -I- a3 , 


r = f /^1 m ^4 _ ^01 \ 

u f, 2/2; ’ 


which to the order of interest also gives 


f = r I 1 - hr ^4 ■>'01 

Jvr Jtt \ y,, 2 ^ 2 ^ 

This result is in accord with that obtained in Ref. [37]. 


L Jn 


(6.16) 


(6.17) 


(6.18) 


C. Renormalization of the ttN coupling constant qa! f-w 

Apart from and which are similar to those in Eqs. (B3) and (B15) (but 

again expressed in terms of renormalized nucleon and pion helds, and pion mass), the other 
interaction Hamiltonian needed is 


Hi^N = / dx ^(2 di6 - dis) + ^ (sZr, + aV ■ 


(6.19) 
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We find that the contributions of the diagrams in Fig. 6 are given by 


bl = cr ■ kTa , 

(10 a - 1) JoifT -kTa, 


b3 = i 
b4 = i 


.• 9 a 


48/3 

r2 


CT ■ k' 


Jtt 
9a 


4/3 


-‘infj L + 


16/3 

10a- 1 


J< 


01 


CT ■ k 


and in terms of renormalized g'’^ and /^ it must be 

. 9 a 


2/’ 


CT ■ k Tq = bl + b2 + b3 + b4 , 


which leads to the following relation valid to order Q 

.2 


9a ^9a 
In In 


9 777^ ^ 

1 H- — (2 die — dig) 




rJi 


= ^(1 + 

/tt 


9a 

ArrfJ g^/ 

“16 ~ “13 ~ 


13 


3// 


/, 


2 

TT J 




3/, 


/, 


r2 


1 - 


2 m: 


r2 


9a 


d 
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( 6 . 20 ) 

( 6 . 21 ) 

( 6 . 22 ) 


(6.23) 


(6.24) 


(6.25) 


where in the second line, in the terms of order Q^, we have replaced gA and In by their renor¬ 
malized values g^^ and /)), which is correct at this order, and have isolated the Goldberger- 
Treiman discrepancy. The above relations are in agreement with Eqs. (102) and (103) of 
Ref. [17]. 


/V ; 


/ / 


b1 b2 b3 b4 

FIG. 6. Diagrams relevant for the renormalization of gA/ In 


Since In has already been renormalized, we can use Eq. (6.25) to independently renor¬ 
malize gA- We hnd up to order 


9 a = 9 a 



2/ 


r2 


dm — 


r2 


3/ 


dl3 -|- 


4 m: 


r2 


9a 


-d 


16 



2m;^ 

9a 



(6.26) 


As a check of this result, in the next subsection we provide a direct renormalization of gA 
by considering the coupling of the axial held to the nucleon. 
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D. Renormalization of the axial conpling constant qa 


The relevant interaction Hamiltonians are and in Eqs. (B20) and (B30), 


and 


H 


(3)/ _ 


ANN - 


9a 


= -/ dxiV' ■ Ai7*75+ (5ZAr^T ■ Ai7*75 + ■ (9^Fij7*75 )iV^ . (6.27) 


We consider a similar set of diagrams as in Fig. 6, but for the incoming pion line replaced 
by the external held. Their contributions are given by 


1 -I 9 a a 

hi = —TafT ■ Aa , 


b2 = 


b3 = 


9a 

9a 

48/2 

9a 


JqI '^a ^ ' Afl, , 


<413 Tq <T ■ Aa , 


d 


b4 = 5Zn + 2 nfJ^diQj TaCr ■ Aa + —Ta (qq ■ cr - q^cr) ■ Aa , 


(6.28) 

(6.29) 

(6.30) 

(6.31) 


and sum up to gA(yTa/ 2 , with the renormalized axial coupling constant (to order Q"^) ob¬ 
tained as 

i — - „ „ aoi — 77777 -'13 H- z — *^16 I 


9 A = 9 A 


2 / 7 -‘ 3 / 7 ^“ ■ * “7 ’ 

and Qai apart from the Goldberger-Treiman discrepancy, is in agreement with Eq. (6.26). It 
is also in agreement with the results, to order reported by Schindler et al. in Ref. 

The term proportional to ^22 quadratic in q contributes to the nucleon axial radius [38]. 


^ / 
> 


d1 


d2 


d3 


d4 


d5 


d6 


FIG. 7. Pion-pole diagrams. 


E. Renormalization of pion-pole contributions 


We examine the pion-pole contributions illustrated in Fig. 7. We obtain 

q - cr 


dl — — — A -a 

~ O -^a 4 2 1 r2 

2 


9a 


Jr 


d2 + d3 = ( -m;2 /4 + ■ q 

9a 


q-o- 

g2 -H ml 

qcr 


:Ta 


= T72 ^ 2 , 

8// m. 


Ta 


(6.33) 

(6.34) 

(6.35) 
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(6.36) 


d5 = 


d6 = 


9a j K q ■ 


48/, 


+ m 


— (2 die ~ di%) + 


dA 

■ 4/2 


2 U + 


16/^ 

10a- 1 


Jf 


01 


Aa ■ q --IT Ta ■ 


-|- m. 


r2 


Their sum reads 


dl + --- + d6 = 


1 - Joi - Ji3 + 


2 f, 


r2 


3/ 


9a 


2 m 


r2 


9a 


d 
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X Aa ■ q ^ T„ 

qZ _|_ 


(6.37) 


(6.38) 


and therefore the renormalized follows exactly as in Eq. (6.26), including the Goldberger- 
Treiman discrepancy. The renormalized (single-nucleon) current is then given by 


: _ 9a , 9a „ q-cr 

35,a — ^ (TTaA- ^ ^ o 2 ^<2- ^ 


2 2 -I- m', 

and this current is conserved in the chiral limit (m^ —)■ 0), since in that limit g\ = gA- 


(6.39) 


F. Renormalization of OPE axial charge 


We begin by discussing the non-pion-pole contributions illustrated in Fig. 8. In diagrams 
g2, g4, g6, g8, gll, and gl4, the solid dot represents the interaction —5m — 4m()^ ci, where 
5m is the nucleon mass counter-term. The contributions associated with diagrams gl-g2, 
g3-g4, g5-g6, g7-g8, g9-gll, and gl2-gl4 represent the renormalization of nucleon external 
lines and, with the choice of 5m in Eq. (6.8), they are seen to vanish. 

Next, the solid square in diagrams gl6, gl8, and g20 represents the interaction 


= - / dx 


dx 


mJ'J' IA ^ 5Z^ 


// ' 2 
(/a -h U) , 


n 


+ 


(n^ ■ n" + 9V’’ ■ 9^77’’) 
5m‘i 


5Z^ 


77 ■ 77 


with vertex (in the convention of Appendix B) 


(0 I I ki, ai; k2, as) = + dz}j (caicas - ki ■ ks) 


2 m'' {l^ + I a) r2x'7 A 2 

-h ml 5 Z^ - 5mt 


P 

J TT 


(6.40) 


(6.41) 


With 5Zt, and 5ml. given in Eq. (6.8), the contributions of diagrams gl5-g20 cancel out. 
The remaining loop contributions in diagrams g21-g29 are given by 


p,g(g21) = p<;‘'(al)^(l - 10a)J„, , (6.42) 
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FIG. 8. Half of the possible time-ordered non-pole corrections to the OPE axial charge at order 
Q. Nucleons, pions, and axial fields are denoted by solid, dashed, and wavy lines, respectively. See 
text for further explanations. 



pU(g22) =pj^„^(al) (1 4a)Joi , 

(6.43) 


pU(g23 + g24) = J(al)2^^2'^i3 , 

(6.44) 


P^;l(g25 + g26) = . 

(6.45) 


pS(g27 + g28 + g29) = . 

(6.46) 

while those in diagrams g30-g32 vanish identically. Here P 5 „^^(al) is dehned as 

( s') f 

Finally, one needs to include the contributions due to the interactions 
and 

Hi^NNA = - (SZn + 6 ZJ2) ^ J d^N^'Ao ■ (r x tt^) yoiV’’ , 

in the OPE axial charge, which simply lead to the correction of order Q 

in Eq. (4.3). 
in Eq. (6.19) 

(6.47) 


2 ^ 1 

2 6Z]\f + 6Z^ + (2 die dig) (al) . 

(6.48) 


9a 


Thus, the sum of the order Q corrections to the axial charge from non-pole contributions, 
denoted as p 5 ^^(npp), reads 





— 


9a 

12 /^ 


Ji3 + 2 6Z^ 


N 


23 



+ H- — (2 die ~ dis) 

9a 

which, which after insertion of 6 Zn and SZ.,^, is expressed as 


(6.49) 


pS(npp) = i 


.• 9 a 


2m;2 

fr2 


Tl X T2)^(T2 ■ k2 — 

2 ml^ 


5 , 5 ^ 7)42 

•^01 — 7rir-7Tai3 




6/. 


r2 


I 4 + 


(2 die — d 


18 J 


(6.50) 


/r 9 A 

where the bare qa and have been replaced by their respective renormalized values—this 
replacement is correct to the order of interest here. The complete non-pole axial charge, 
denoted as p 5 ™(npp) below, results from the sum of the leading-order contribution in 
Eq. (4.3) with the ratio Qa/ fI replaced by its renormalized value 

1 


9a 


f2 fr2 
J 7T J TT 


2fr2 3J.2 

and the contribution p 5 ^](npp). We obtain 


doi + o H- U -(2die — dig) 


f: 


9a 


(6.51) 


P?™(npp) = (ti X T2 )„ 0-2 ■ k2;^ ( 1 + 


. 9 a 


81. 


I )2 

UJ 2 


8 f:, 


r 2 


T 9 A J 
9oi - -'IS 


2/, 


(6.52) 


The diagrams describing the pole corrections are illustrated in Fig. 9 (only representative 
diagrams for each of the relevant classes are drawn for brevity), and are similar to those in 
Fig. 8. A slightly more complicated analysis along the lines illustrated above leads to a pole 
OPE axial charge, denoted p 5 ^](pp), given by 


PSr(PP) = ('^1 X 'r2)a 0-2 ■ k2;^ ( 1 - 


. 9 a 


8 f:. 


, )2 

UJ 2 


81 . 


r2 


T 9 A J 

901 - 9i3 


2/. 


(6.53) 


The sum of the npp and pp contributions evaluated in dimensional regularization is 


P?™(npp + pp) = (ti X T2)„ 0-2 ■ ks;^ ( 1 - ^ J, 


, 9 a 


8 f:^ 


/ )2 
U 2 

3 -771^2 
''O r2 




fr 
J TT 




\^2 fr 2 - 


9^ 4 - X 


(6.54) 


There are additional loop corrections to the OPE axial charge, see Fig. 10. Their contribu¬ 
tions are obtained as 


T" 2 

pS(fl + f2) = , 

40 


p^:^(f3+f4) = 
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r2 
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(6.55) 
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where P 5 ^^^(al) is again defined as in Eq. (4.3), except that gA and are replaced by their 
renormalized values and /)). The loop function has been defined in Eq. (4.7), 

while I^‘^\k) and L{k) read 

dp 


I^^\k) = 


(27r) 


P /(x;-,x;+) , 


(6.57) 


24 




h11 h12 h13 



FIG. 9. Representative diagrams for each of the relevant classes contributing to pole corrections to 
the OPE axial charge at order Q. Nucleons, pions, and axial fields are denoted by solid, dashed, 
and wavy lines, respectively. More than a single time ordering is shown for some of the diagrams. 
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FIG. 10. Additional loop and tree-level corrections of order Q to the OPE axial charge. Nucleons, 
pions, and axial fields are denoted by solid, dashed, and wavy lines, respectively. Only a single 
time ordering is shown for each topology. See text for further explanations. 
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Evaluation in dimensional regularization leads to 
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with as given in Eq. (4.19). We also need to account for tree-level contributions of order 
Q originating from the vertices 2t[NN and NNttA^ in Eqs. (B14) and (B29), denoted by 
the solid diamonds in Eig. 10. They can be written as 


p^^](f5 + f6) = 2 (di kl + d2kl + ds + d^ 


+ *^5^2,0 CTi ■ (q X k 2 ) <72 ■ k 2 ^ , (6.61) 

^ Jtt ^2 

where we have introduced the following combinations of LECs 

di = 2 d2 + de , (6.62) 

(^2 = 4 di -|- 2 (^2 + 4 ds — de , (6.63) 

d3 = -2d2 + dQ, (6.64) 

d^ = 4 di -|- 4 d2 “1“ 4 d^ -|- 8 d^ , (6.65) 

ds = di5 -|- 2 d23 . (6.66) 


The divergent parts of the dj’s (and hence dj’s) have been identihed in the heavy-baryon 
formalism, without considering any specihc process, with the background-held and heat- 
kernel methods, see Ref. [39] and references therein. We report below the expressions for 
these divergent parts from Table 4 of that work: 


di — -I- d\{p) , 

J TT 

where, in the conventions adopted in the present work. 
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(6.71) 


and 05 is from Eq. (B13) of Ref. [39] which corresponds to our choice of operator basis in 
For the combinations di above we obtain 
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and ds = dj and d^ = d^. We observe that the divergence proportional to from loop 
corrections in pPjj^^(npp + pp) cancels exactly that present in fl + f2. Next, the divergent 
part of di cancels exactly the term proportional to kf d^ present in fl + f2 and f3 + f4. The 
divergent parts of d 2 and d/^ are the same, and therefore can be reabsorbed in the LEG Z 2 
mnltiplying the contact term 02- Those of d^ and d^ vanish, which is consistent with the 
fact that there are no divergencies proportional to or in the operator mnltiplying d^. 

Combining Eqs. (6.52), (6.53), (6.59), (6.60), and (6.61), we then End that the renormal¬ 
ized OPE contribntions np to order Q inclnded read as 
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(6.75) 


G. OPE axial current 


The loop corrections to the OPE axial cnrrent are shown in Figs. 9 and 11. Those 
associated with panels hl-hl7 are easily seen to vanish, while the contribntions of diagrams 
ml-m2 are obtained as 


5 

jS(ml) = Ju [9r2,ak2 - (Ti X T2)„ (<Ti X k2)] 0-2 ■ k2 ^ , 

jS(m2) = - q-j^ll(ml) , 

In dimensional regnlarization we find the finite result 
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No renormalization is necessary in this case, since loop corrections to diagrams dl-d2 of 
Fig. 4 enter at order Q", and are beyond the scope of the present work. 


VII. DISCUSSION 

In this section we report the complete (and renormalized) expressions for the weak axial 
charge and current operators, compare these expressions to those obtained by the authors of 
Ref. [18], and discuss current conservation in the chiral limit. For simplicity, the superscript 
r has been removed from the pion and nucleon masses m^r and m, the nucleon axial coupling 
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FIG. 11. The only non-vanishing loop corrections to the OPE axial current. Nucleons, pions, and 
axial fields are denoted by solid, dashed, and wavy lines, respectively. Only a single time ordering 
is shown for each topology. 


constant qa-, and pion decay constant /^r. However, all these quantities are understood to 
have been renormalized. 

The one-body operators and two-body contact operators are those listed, respectively, in 
Eqs. (3.13)-(3.14) and Eqs. (5.1)-(5.2), and in Eqs. (4.20) and (5.4), while the two-body 
operators involving OPE, TPE or MPE, and short-range terms follow in the next subsection. 
Relativistic corrections (proportional to 1/m^) in the one-body axial charge are neglected, 
those in the one-body axial current (proportional to 1/m^) are retained in Eqs. (5.1)~(5.2), 
since they are known to be important in weak transitions such as the proton weak capture 
on ^He at low energies [33]. 


A. Two-body axial charge and current operators up to one loop: summary 

The (renormalized) OPE contributions to the axial charge are given in Eq. (6.75), while 
those corresponding to the axial current read 
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The TPE axial charge, and MPE and short-range axial current can be written, respectively, 
as 
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with with Sj dehned as in Eq. (4.16), and 
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and the loop functions Zi and Wi are listed in Appendix D. 


(7.4) 


(7.5) 


B. Current conservation in the chiral limit 


In the chiral limit (m^r —)■ 0) the axial current is conserved and 

q- j5,a =[H , P5,a] , (7.6) 

with the two-nucleon Hamiltonian given by 

^ = + , (7.7) 

where the superscripts denote the power counting, the are the two-nucleon potentials 
dehned in Sec. Ill, and the kinetic energy (in momentum space) is 


b = (27r)^(5(p2 - P 2 ) + (1^2). 

2m 


(7.8) 


Here, the potentials and axial charge and current operators (including the axial coupling 
and pion decay constants and LECs entering them) are to be understood in the chiral limit. 
Order by order in the power counting, current conservation implies the following set of 
relations 
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where we have only kept up to terms of order Q^. Note that the commutators implicitly 
bring in factors of Q^. The hrst of these relations is obviously satished, see Eqs. (3.14) 
or (6.39). The second relation has 

q ■ = -^ ri,a ki ■ Ki 0-1 ■ Ki + (1 ^ 2) , (7.13) 

where is given by the sum of the contributions in Eqs. (5.1) and (5.2), and is also 
satished. The left-hand-side of the third relation has 
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and this matches the hrst commutator on the right-hand side, T*- with given 

by 
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i.e., the sum of terms al and a2 in Eqs. (4.3) and (4.4). There are additional contributions 
to which arise from non-static corrections to the denominators involving pion energies 
in the diagrams illustrated in Fig. 12, where the crossed circle (cross) means that the these 
denominators are expanded as indicated in Eq. (3.3) to order Q (Q^) beyond the leading-order 

static term. These contributions are needed in order to satisfy the commutator 
but have been neglected in the present work. 
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FIG. 12. Illustration of some of the non-static corrections to the axial current ignored in this work. 
Nucleons, pions, and axial fields are denoted by solid, dashed, and wavy lines, respectively. See 
text for further explanations. 


Lastly, we consider the fourth relation, Eq. (7.12). The axial current obtained here 


is in the static limit, and one expects q ■ jg^^ to satisfy the commutator 
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where the loop functions R'^'^\k) and R\j {k) in the chiral limit read 
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and the ... indicate a term proportional to R kj, which vanishes when inserted in Eq. (7.16). 
The current-conservation constraint is seen to be satished by noting the only non-vanishing 
contributions to q ■ are those due to diagrams e4, e5, elO, el5, el6, and el7 in Fig. 4, 
proportional to the combination of coupling constants g\/f^. In particular the contributions 
of the purely irreducible diagrams e4, e5, el5, el6, and el7 combine to give 
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with the remaining “missing” term being provided by q-j5^^(el0). The remaining commuta¬ 


tor 




has a factor 1/m, and therefore non-static corrections need to be included 

■ J 

in if the latter is to satisfy the complete Eq. (7.12). These corrections have been ignored 
in the present work. 


C. Comparison 

We compare the one- and two-body axial charge and current operators derived here with 
those obtained by Park et al. in Refs. [18] and [33] in the heavy-baryon (HB) formulation of 
covariant perturbation theory. The one-body axial charge and current operators at leading 
order in Eqs. (3.13) and (3.14) are the same as those listed in Eqs. (Bl) and (A3) of Ref. [33], 

/ Q \ 

except for the pion-pole contribution to jg ^ q which, while nominally of the same order (Q”^) 
as the non-pole contribution, is nevertheless suppressed at low momentum transfer q and 
is therefore ignored in Ref. [33] (we note incidentally that in that work ki = —q, i.e., the 
opposite convention adopted here). Of course, this pion-pole contribution is crucial for 
current conservation in the chiral limit. We have neglected the 1/m^ relativistic corrections 
to the leading order axial charge. They are retained in Eq. (17) of Ref. [33]. However, the 
1/m^ corrections to the leading order axial current in Eq. (5.1) are in agreement with those 
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given in Eq. (A3) of Ref. [33], except for the last term proportional to q (cti ■ q), which is 
again ignored in that work. 


FIG. 13. Feynman amplitudes contributing to the one-body axial charge at leading order. 


Before moving on to the two-body contributions, it is worthwhile discussing how the 
one-body axial charge operator emerges in covariant perturbation theory. The relevant 
interaction Hamiltonian densities are 
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(7.22) 


= N{x)t ■ di7r{x)x''y^N{x) , 

(7.23) 


where all helds are in interaction picture. The R-matrix elements associated with the Feyn¬ 
man amplitudes in Fig. 13 are given by 
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where 7 = a or b, T denotes the usual chronological product, and |p. A) and |p'. A') are the 
initial and hnal nucleon states with momenta p and p' in spin-isospin states xa and xx'^ 
respectively. Then for 7 = a we obtain 
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where we have considered the leading order in the non-relativistic expansion of the nucleon 
matrix element. Since in the interaction picture the conjugate held momentum nc(a;) = 
d^nc{x), it is easily seen that (see also Ref. [40]) 
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with the Feynman propagator dehned by 
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The T-matrix element Tfi obtained from Sfi = —i {2n)‘^ 6{p' — p — q)Tfi reads 


Tff = cr ■ (P + P) TcXa ( 1 + 




ml + q^ — Qq — ie 


(7.28) 


where the term proportional to go = Po ~ Po is suppressed by in the power counting. The 
leading order term leads to the axial charge operator in Eq. (3.13). A similar analysis shows 
that the leading-order contribution to vanishes. 

As already noted, the interaction Hamiltonian in Eq. (2.11) contains no direct coupling 
of A° to the nucleon. However, diagrams of the type illustrated in Fig. 13 are not considered 
in Refs. [18, 33]. It would appear that their contribution is accounted for by retaining the 
term —i6cdS{x — y) in Eq. (7.26), which effectively leads to a direct coupling between A° 
and the nucleon. 

Turning to the OPE contributions at tree level, we hnd that the order Q~^ contribution 
to the axial charge, p^^a \ Eq. (6.75) reproduces the corresponding contribution, given by 
Eqs. (B2), (B3), and (B5) of Ref. [33] with FY{t) = 1, while the order contribution to 
the axial current, in Eq. (7.1) is the same as in Eq. (A5) of Ref. [33]. We stress again 
that, while diagram a2 in Fig. 2 is not explicitly considered in Refs. [18, 33], the OPE axial 
charge operator derived there has the correct strength. The contact terms contributing to 
the axial current in Eq. (A6) of Ref. [33] can be reduced through Fierz identities to the 
form given in Eq. (5.4). 



n 1 




n3 


FIG. 14. Diagrams contributing to the axial charge (nl-n2) and current (n3) at order Q consid¬ 
ered in Ref. [18]. Nucleons, pions, and axial fields are denoted by solid, dashed, and wavy lines, 
respectively. Only a single time ordering is shown for each of the possible 12 (nl) and 60 (n2 and 
n3) cross-box topologies. 


Next we consider loop corrections to the axial charge. The contributions of c3-c4, c7-c8, 
and c9-cl2 in Fig. 2 are found to vanish in both approaches, here and in Refs. [18, 33]. The 
contributions of diagrams cl and c2 are the same as those for A*'°^(a -|- b) in Eq. (93) of 
Ref. [18]. The contributions of diagrams c5 and c6 are different from those for A^^\c + d) 
reported in Eq. (94) of Ref. [18] because of the different treatment of reducible topologies 
for these types of terms. Indeed, if only the (irreducible) cross-box topologies are retained 
for diagrams c5 and c6, as illustrated in Fig. 14, then the resulting operator is the same as 
in Eq. (94). The OPE axial charge operator in Eqs. (74) and (90) of Ref. [18] reads in our 
notation 
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Provided we define 


d; + d; - dl - = -(ITgi + 4)/(144»VJ) - c; , 

the expression above is in agreement with our Eq. (6.75) in the limit q = 0 (or ki = —k 2 ) 
which is assumed in Refs. [18, 33], except for the term proportional to rri^ in the first line. 

Lastly, the term proportional to the LEG C 3 in Ref. [18] (in the HB formulation) is given 
by 

[Z}/3, [D^,Dp]]N, 

J TT 

which can be re-expressed as 

,^N[D>’,F„*,]N + ... . 

and matches the term proportional to dg in fhe HB limit of [27] —in the relation above 

n" is the velocity, n" = ( 1 , 0 ). 

Moving on to the loop corrections to the axial current, the sum of the contributions 
due to diagram ml of Fig. 11 and diagram el5 of Fig. 4 gives the same expression as in 
Eq. (A7) of Ref. [33], provided the parameter a in the 37 rH vertex of diagram el5 is set 
to 1/6—the authors of Refs. [18, 33] use the exponential parametrization for the pion field. 
The irreducible contributions of diagrams el and e4 in Fig. 4 are the same as reported 
for, respectively, A“2(27r:b) and A“2(27r:a) of Eq. (A13) of Ref. [33], while the contributions 
associated with the cross-box topologies of diagram e 8 in Fig. 4 and illustrated in panel n3 of 
Fig. 14, lead to the expression for A“2(27r:c) in Eq. (A13). Non-vanishing pion-pole diagrams 
e2, e5, e9, elO, el 6 , and el7 as well as diagrams e20-e21 (e22 and e23 vanish) in Fig. 4 have 
not been considered in Refs. [18, 33]. In particular, because of this incomplete treatment, 
loop corrections to the axial current are a-dependent in Refs. [18, 33]. Furthermore, the 
current is not conserved in the chiral limit. 

Finally, the OPE axial current at tree-level listed in the recent Ref. [41] (and including 
pion-pole contributions) is different from that obtained in the present work in Eqs. (5.5)- 
(5.6). Moreover, it is not conserved in the chiral limit. 


VIII. CONCLUSIONS 

In the present work we have carried out an analysis of the weak axial charge and current 
operators in a two-nucleon system up to one loop (i.e., including corrections up to order Q 
in the power counting) in yEFT. The formalism used in the derivation is based on standard 
TOPT, but accounts for cancellations between the contributions of irreducible diagrams 
and the contributions due to non-static corrections from energy denominators of reducible 
diagrams. A detailed comparison between the results of this work and those of the early 
studies of Park et al. [18, 33] in the HB formulation of yEFT indicates that there are 
differences in some of the loop corrections and in the renormalization of the OPE axial 
charge, the former due to a different prescription adopted by the authors of those papers, 
one in which only a subset of the irreducible contributions are retained in the perturbative 
expansion—for example, in the case of box diagrams, only cross-box ones are considered. 
Furthermore, while the contribution illustrated by panel el5 in Fig. 4 is accounted for in 
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Refs. [18, 33], additional ones involving three- and four-pion vertices, such as those in panels 
e5, el6, and el7, have been ignored. As a consequence, the one-loop axial current derived 
there depends on the parametrization of the pion held—it is a-dependent—and, furthermore, 
is not conserved in the chiral limit. 

The order Q loop corrections in the axial current are hnite, consistently with the fact 
that there are no contact terms at this order. There is a single LEG (denoted as Zq here 
and as dn in Ref. [33]) which enters at lower order Q^. On the other hand, four independent 
LECs (denoted as Zi, with i = 1,... ,4) multiply contact terms in the axial charge at order 
Q, two of which are needed to reabsorb the divergencies from loop corrections in the TPE 
axial charge. The loop corrections to the OPE axial charge instead lead to renormalization 

•—' /Q\ 

of di which is expressed as linear combinations of the LECs di in the Lagrangian—some 
of these di having been determined in hts to nN scattering data [42]. The LEG Zq has been 
recently fixed by reproducing the empirical value of the Gamow-Teller matrix element in 
/d-decay [36]. However, that calculation ignored MPE loop corrections in js^, and therefore a 
rehtting of zq will be necessary. Most calculations of nuclear axial current matrix elements, 
such as those reported in Refs. [33, 43] for the pp and p^He weak fusions and in Ref. [36] 
for muon capture on and ^He, have used axial current operators up to order (one 
exception is Ref. [44], which included effective one-body reductions, for use in a shell-model 
study, of the TPE corrections to the axial current derived in Ref. [33]). Lastly, there remains 
the problem of determining the Zi^s in the contact axial charge. It should be possible to £x at 
least some of these LEGs by studying muon capture in the few-nucleon systems, for example, 
by reproducing data on angular correlation parameters for the process ^He(/i“, z/^)^H [45], 
or cross sections for transitions from the bound state to breakup channels, such as the ^H-n 
two-body breakup, for which data are available [46]. 

On a longer time scale, it should be possible to use the weak axial operators constructed 
here in quantum Monte Garlo (QMG) [47] calculations of /3-decays and electron- and muon- 
captures in heavier nuclei with mass number A > 4 (see Ref. [48] for an earlier study 
of these processes in ®He and ^Be in the conventional meson-exchange framework) and of 
neutrino inclusive cross sections off light nuclei at low energy and momentum transfers [49]. 
As a matter of fact, the very recent development of “realistic” and mildly non-local chiral 
potentials in conhguration space [50], in which QMG methods are presently formulated, 
makes it possible to carry out these calculations in a consistent yEFT framework (i.e., 
chiral potentials and currents), and hence offers the opportunity to provide first-principles 
(and numerically exact) predictions, rooted in QGD, for the rates and cross sections of these 
weak processes. 
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Appendix A: Chiral Lagrangians 


We adopt the notation and conventions of Ref. [27] for the various helds and covariant 
derivatives, which we summarize below: 


U 

u 


1 + ■ 77 - 

Jtv 

Vu = i + 


1 % cx ry 

- TT^ — - TT^ T • TT -I -TT^ H- 

2 p ^ ^ “T . . . , 


8 /^ 


2 /, 


-T ■ 77 — 


8 /i 


2 7 (8 a - 1) 2 , (32 a-5) 4 

77- „„ 77 T ■ 77 H-—; 77 + . . 


16 


128/4 


(Al) 

,(A2) 


Uf, = i[u\df,-irf,)u-u{df,-ilf,)u'^] , (A3) 

D^U = d^U-ir^U + iUl^ , (A4) 

D^N = {d^ + T^)N = df,N + ^ [u^id/, - m + m {d^ u^] N , (A5) 

Fitu = , (A 6 ) 

-i[r^,r^] , r^ = v^ + a^ , (A7) 

A/jy dfxlu d^il^ i [ /^, /i/ ] , VOj^ , (A 8 ) 

X± = Xu±ux^ u = ml. (U^ ±U) . (A9) 


The parameter a is arbitrary because of the freedom in the choice of pion held—the only 
constraint is that U be unitary with det 17 = 1. Common choices are a = 0 and a = 1/6 
corresponding, respectively, to the non-linear sigma model U = {a + i t ■ tv) / with a = 
— 77^ and to the exponential parametrization U = exp(iT ■ 77 //.^-). In the following 
we consider only the coupling to the axial-vector held; further, we ignore isospin-symmetry- 
breaking ehects as well as the coupling to the isoscalar component of the axial-vector held, 
and hence 

r^, = , (AlO) 

F^u = ■ ('9/.A^ - dyA^ + X A,,) , (All) 

F!^v = ' ("^mA,. - '9^A^ - A^ X A^) . (A12) 

Inserting the expansions for U and u and keeping terms linear in the axial-vector held, we 
hnd: 




D^U 

D^N 


1 f ^ a A . da — 1 ^ 

+ ■ A;. + ^ [(r X 77 ) X 77 ] ■ A^ + ... , 


l T 


J 7T 


^ 1 A 

2 /r ^ 


1 0,1 A , 

^ 77 ■ d^TT + — TT ■ Af, + . . 

Jtt If 




77 


+l 


_.( 8 a - 1)_2 


16 /^ 


77 duTT ■ (77 X r) -I- 


1 “ "72 ) ('^ ^ '^) ■ A^ 

J TT 

N , 


(A13) 

(AM) 


(A15) 
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F;1u = ^ (t X tt) ■ 




F- = 


2 /^ 


T X 77 X 77 


F 


IIU 


77 


X+ = (^2 -—)+.. . , 

2i 2 

X- = ■ 77 

In 


(A16) 

(A17) 

(A18) 

(A19) 


where F^^, = dfj_Aiy — di^A^ and the ... denote higher powers of the pion held than shown. 


1. 77A sector 


The nN Lagrangians up to order read: 

4‘t = , (A20) 

£S = Z<;.]VOf'iV, (A21) 

i=l 

23 

= (A22) 

^=1 

with the operators and dehned as in Ref. [27]. Here gA is the nucleon axial 
coupling constant, and the q and di are LECs. Below, the 7 ^, 75 , and are 7 matrices 
and combinations of 7 matrices in standard notation [51], and is the Levi-Civita tensor 
with = + 1 . 

In terms of the expansions above, is given by 


^■jtn — 


m 


4/^ 


T X 77 ■ WTr 


9A 


+ J^{4a- 1) T ■ 7777 ■ (^7775 + 4A—^pi7r^^77 ■ (77 X r) 

4 /i? 16 n 

1 


2f-, 

(1 - 8 a) 


1 -TT^ 

P 

J TT 


T -077 


2 /. V /, 


1 “ z w 

1 - ^77 fr X 77 


+ ^ T ■ A 75 + [(t X 77 ) X 77 ] ■ 

2 4/,^ 


■75 


N 


(A23) 

where = 7 ^( 9 ^ and = 7 ^A^. The operators 0^1 the Lagrangian are expressed as 
(below the notation x+ = X+ ~ (x+)/2 is used, where (...) implies a trace in isospin space) 

of = (x+> ^ 4 mj ([1 - |7]) , (A24) 

of = + h.c. —> ■ 9oJr - ■ Ao 

0 777 Jtt 

H—V f “ '^ 0 '^ ■ "y” 7 V* (A25) 

777/^ VA / 
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C'? = ^ (9^77 ■ (9^77 ■ , 

^ In Jtt 


of'’ = ^[u^,, u^] (y^'" 


of'’ = x+ 


0 , 


(^( 2 ) _ 
*^6 “ 


of’ = 


- F+ - 

8 m 

8 m ^ 

’(3) 


4m/7r 
-)■ 0 , 


f2 ' p 

J 77 J 

(r X 77 ) ■ (9^A^, , 


1 2 

^ '^ ■ ( “72 ^ 1 

1 


(A26) 

(A27) 

(A28) 

(A29) 

(A30) 


while those in the £^j(r Lagrangian reduce to 


O® = 


, [D ^, + h.c. —> Y'^' {~ ^ <9o«9^77 + A^ X (9o(9^77 


-(9oA^ X (9^77^7° , (A31) 

of = > [O^ > + h.c. —> y'^ ' {~ ^ d^doTT + A^ X dod^Tz 


( 9 ( 3 ) ^^ 

^ 12 m3' 


fn V /tt 

-(9^Ao X (9^77^7° , 

[u^ , , Up\]D^^''P + h.c. — > Y'^ ' i~ ^ *^0^ + Ao X (9, 

I TT ' / 77 


77 


-(9oAo X (90 77^7° , 


= 


2 m 


Yf,UuUa)Dp + h.c. —)■ 0 , 


of = 


2 m 


[X- , Uf,]D>^ + h.c. 


4m^ 

'If 


TT X ( —(9o 77 - Ao 

J 77 


7 


( 9 ( 3 ) _ * 

IXg — 


2 m 

(3) ^ 


= 


2 m 


P',F+l£>‘' + h.c, ^S'7„7“, 

P', {F+)]D-' + h.c. ^ (I , 


q(^) — ^ 


2m 


{Ffua)Dp + h.c. —^ 0 


of = — 

2 m 

<^10 = \ 7^75 {u-u)Uf, —^ 0 , 

C’n = \ 7^75 {u^f u'' —^ 0 , 


/ira/3 


{Ff) Dp + h.c. —> 0 


( 9 ( 3 ) ^_ _ 

8 m 2 

( 9 ( 3 ) __ 

“ 8 m 2 


7^75 {uxuj) Uf, D^'' + h.c. —)■ 0 


fxb (uYu) U\ D^'' + h.c. —)■ 0 


( 9 ( 3 ) _ 

<^14 — 


4m 


a^''{[Dx , u^] Uy) + h.c. 


1 

Iff: 

(9oAi ■ djifj o-*^7° , 


— f —(9o(9j77 ■ (9,-77 — (9o(9j77 ■ A,- 


(A32) 

(A33) 

(A34) 

(A35) 

(A36) 

(A37) 

(A38) 

(A39) 

(A40) 

(A41) 

(A42) 

(A43) 

(A44) 
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oSJ = ^ , ux\)D^ + h.c. —)■ y ~ 

-Ai ■ dodjiT^ o-*^ 7 ° , 

<^16 = \ 7^75 (x+) “m —> 2 ml T ■ (^ - ^(9*77 + Aj^7*75 , 

Ofi = ^ 7^75 (X+ u^) —^ 0 , 

<^18 = \ 7^75 [Df, , X-] —^ ^ ’ 

oS? = ^ 7^75 > (X-)] —^ 0 , 

oS = -^ 7^75 [F+ , «a] + h.c. 0 , 

Og^ = ^7"75[F+,«1^0, 

Og^ = ^ 7^75 [D'^ , -^\t- 7*75 , 




■-e*"“^9,7r-F„;3 7*75 
Jtt 


(A45) 

(A46) 

(A47) 

(A48) 

(A49) 

(A50) 

(A51) 

(A52) 

(A53) 


Several comments are now in order. First, the expressions above for and £gg 

retain all terms relevant in the present study. Typically, these include at most three pion, 
two pion, and one pion fields for n = 1, 2, 3 in respectively. In some instances, for 

/qA 

example in , terms with two pion fields are also considered for reasons having to do with 
the treatment of tadpole-type contributions (see below). The Lagrangian 
conveniently be expressed as given in Eq. (2.4) with the quantities r°(n), A)j(n), and A(n), 
dehned in Eqs. (2.5)-(2.6), given at leading order by 


^ ^ ^)a7°, 

4 11 16 

^a(O) = fl- T-a 7*75 + ^ (4 a -1 )(t- 77)77^7*75 , 


^(1) = Y 9'*75 + ^ ^ 

+ /4 [(t X 77) X 77] ■ Ai7*75 ; 


(A54) 

(A55) 


(A56) 


at next-to-leading order by 


ra(l) = (^1 - + ^(4a-l)(T-77)7ra7°75-2 ^^^^^ A° , (A57) 


^ 7^)ay+ 5'9Va-2^A),- ^ (r X Aj)^a*^ 
^ In Jtt Jf Jtt 

+ ^ ('^ X dj7r)^a^^ + (tt x T)a7* , 


(A58) 
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^(2) = yT ■ Ao 7°7^ ^ ^ ~ ('^ X ■ Ai 7* + [(t x tt) x tt] ■ Ao7°75 


+4m^ c, ^ j + ^ (r X TT) . 5.A, a- ; 

and at next-to-next-to-leading order by 


(A59) 


- 0 /o^ _ ^2 / ^ /i . ^4 


r::( 2 ) = 


mfn \U J Jtt 


— (t X diTz)^ -{rxAi 
J-K 


a 


Oi 


+2 


+2 


di -\- d2 

J TT 




X (9*(9i7r) -I- (r X (9V) 6 


7 


di -|- (i 2 + ^3 


(r X tt)^ - ^ (r X <9Vi7r)^ + (r x diA^)^ 


f. 


mi , , 0 (ii 4 - di 5 


-4 4^ (t X 7r)^7 + 




A 


/, 

y (9i7ra (T*^ (9iAj'a CT*^ 

In 




+Aj^a 0 '*'^ <9 


(^2 


+ , 

Jtt 


r 


(A60) 


Al(2) = — ^ + |!- (r X A„)„^“ -2^{-rx d,A')J^ 

m Jtt Jtt Jn 

-27(^xs-A„)„/-7(^xF»y7 

J TT J 7T 

-h^(9oAj-aO-*^7° - ^(9jAo,aO-*^7° 


/. 


m: 


/. 


^ (2 di 6 - dig) Ta 7*75 + 2 ^ e*^'=°Ffco,a7,75 , 

Jn Jn 


(A61) 


2^(3)= , (r X tt) ■ ((9oAi - diAo) o-°* _ 2 _ ^ + 4 


4m/„ 




(r X Ao) ■ ((9*(9^77 -h m^7r) 7 “ 


2 7 

-h4 dg ^ T ■ (tt X Ao) 7 ° -h y (r X tt) ■ (9Tj0 7° + 2 die r ■ Aj 7*75 
Jn Jn 

d^‘2 0(2 T-' ^ 

-F—T ■ (9 Fi^7 75 . 


(A62) 


Second, the various derivatives act only on the held to their immediate right, for example 

^^ ^ ^ ^ ^ ^ _ 

doTz ■ Ao means (( 9 o 7 r) ■ Ao- However, the symbols d i = d i — d i and d i = d i + d i m. 
Eqs. (A25) and (A60)-(A61) denote derivatives acting only on the right and left nucleon 
helds, respectively. 

Third, the power counting Q” of counts powers of derivatives of the pion held (or of 
pion mass factors) and factors of and its derivatives (note that is counted as being of 
order Q). However, the Lorentz structure of the terms may lead to additional suppression. 
For example, in a term like 


4/) 


-{t X tt) ■ (9o7r7° , 
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is of order Q, but a term like 


9A 

2/. 



T 


'9o7r7°75 , 


which is nominally of order Q, is in fact of order since couples the lower to the 

upper components of the spinors, and therefore involves the three-momenta of the initial 
and final nucleons (of order Q). We have taken advantage of this suppression in some of 
the terms in by retaining only the diagonal piece in their Lorentz structure, for 
example in term 

Fourth, time derivatives of the nucleon fields in and are removed by making 
use of the equation of motion (to order Q) 

doN = -imdN + i7°r°(0) doiia -f f 7°Aa(0) diHa -f * 7 ° A(l)] N , (A63) 


implying that 


dlN=—w? N — imd\ - ■ - — im [•••] 7 °A^ 

mgA 


(A64) 


=—m N -I- 


TTt 

T ■ diTT 7 * 7 ^ + mgAT ■ Ki 7*75 - —r ■ ( Aq x 77 ) 7 ° N , (A65) 
/tt Jtt 


where in the second line we have ignored non-linear terms in the pion field, since they do 
not contribnte to the order of interest here. 

f 3) 

Fifth, donble time derivatives of the pion fields in are removed by making use of 
the equation of motion, see Eq. (A72) below. Terms containing both one time derivative 
and one space derivative of the pion fields have been rewritten by integrating by parts. For 
example, in d^jj- a term like 


^di+^ A (r X (9o(9*7r) ■ (9*77 A , 

J TT 


can be re-expressed, modnlo a total divergence, as 
di + d2^ 


- 2 - 


f2 
J TT 


-N 




(r X (9 o77) ■ (9*(9*77 4- (r X don) ■ (9V d 


A 


2. 7777 Sector 


The 7777 Lagrangians np to order read [28]: 

f { D,U (£»■[/)* + x+ ) 


(A 66 ) 


(D>‘U)^) {DM (DM)') + ‘j{DM (DM)') {DM (DM)') + F(x+ )= 


y 
16 L 


2 ( D,U (D^U)^ ){^^) + 2{dudU + xU^xU^)-{X-f-^{ dx) 


+k {{Fd U Fd U^) - ^ (F^Fd + F^FI^^}') 
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+* |< -C (D-U)' + (D^U)' D-U) - 

((x+>" + (X-)" - 2 (xC/t X + Ux^Ux')) 

-2 h,{F;:, Fr + n") . (A67) 

where in the absence of isospin symmetry breaking (which is assumed throughout the present 
work) X is proportional to the identity matrix, namely y = and (x-) vanishes. Further¬ 
more, the terms proportional to the LECs /i, / 2 , I 5 , Iq, and hi do not contribute to the order 
of interest. The symmetric matrices Gab, Gab, Hab, and Fab in the Lagrangian of Eq. (2.4) 
are obtained as 


Gab — 
Gab = 

Hab = 

Fab = 


2 a 


mt 


1 -+ 2 hab ~ 


G 


^C2 + C3 — 


4 q; — 1 


J 77 




ab 


P 

J 77 


NN6ab , 


1 - 

1 - 


8 a 


ml 


2 Q; -|- 1 2 


77 + 2 (I 3 + I 4 ) 


^ab 


mt 


2 /^ 


77 + 2 I 4 ) 6 ab — 
J TT 


2 a- 1 


P 

J TT 


■ T^a^b 


(A 68 ) 

(A69) 

(A70) 

(ATI) 


By retaining only terms linear in the pion held and external axial held, the equation of 
motion implied by is 


d^TT = - [d'-di -h ml) 77 + /^(9oA° -h f^diA". 


(A72) 


Appendix B: Interaction vertices 


In this appendix we report expressions for the vertices corresponding to the interaction 
terms in the Hamiltonian of Eq. (2.11), which we write as 

3 

TFT _ \ ^ r f ZjG) _1_ ZjG) I ZjG) I 'l I TjG) _i_ IjG) _i_ ZjG) i 'i 

“T ^2ttNN + ^3ttNN " J ~r “r ^nNNA “r ^2nNNA ' J 

n=l 

2 


+ E [ (^S”" + n2r'+ ■■■) + (AT + A'”’ + • • •) 


(Bl) 


m=l 


where the superscript n denotes the power counting and the subscript specihes the num¬ 
ber of pion, nucleon, and axial helds entering a given interaction term. We use the following 
notation: A = p ur (A' = p' Pr') are the momentum and spin and isospin projections of 
the initial (hnal) nucleon; ki,k 2 ,... and 01 , 02 ,... are the momenta and isospin projec¬ 
tions of pions 1 , 2 ,... with energies coi, UJ 2 , ■ ■ ■, where Ui = x/kf + mp, q and o denote the 
momentum and isospin projection of the external axial held with energy Uq and its spatial 
and time derivatives expressed as —> i^.^a OqA^^ —)■ —iojqA'^. We also dehne 

P = (p' -|- p)/2 and the (inhnite) constants 


Jn 


dl 

( 277)3 


(B2) 
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1. ttNN vertices 


The interaction terms read 
^(1) 


^ / dxiVr ■ diTvYj^N , 


H. 


2fn 

(2) fi'A 


ttNN 


= ^ / dx Nt ■ 117 ^ 7 '"iV , 
^Jtt J 


H. 


(3) 


m. 


ttNN 


= -j^{2die - dis) J dxiVr ■ (9*777, 


from which the following vertices for pion absorption are obtained 

I hI^nn I k, a) = i-^Ta cr ■ k , 

^Jn 

I I k, a) = -i Ta U(T P , 


(-^'1 l-^;k,a) = i^(2di6 - di8)r„ cr ■ k + 


2 m 
2 


-cr fp - pj 


Sm^f, 

(p' - p) ■ k 2 


(B3) 

(B4) 

(B5) 

(B6) 

(B7) 


2(TPkP 
2 cr ■ k — (k ■ (p' — p) X P 


,(B8) 


where on the r.h.s. of the above equations the 1/\/2 ca normalization factor from the pion 
held expansion in normal modes, the initial and hnal spin-isospin states of the nucleons, 
and the three-momentum conserving (5-function (27r)^(5(p' — p — k) have been dropped for 
simplicity. We will continue to do so in the equations to follow. Vertices in which the 
pion is in the hnal state (pion emission) are obtained from those above by the replacements 
ca, k —)■ —ca, — k. Lastly, only the leading order is retained in the non-relativistic expansion 
of the Lorentz structures associated with the various interaction terms (here and to follow) 
unless otherwise noted. 


2. 2ttNN vertices 


H. 


The interaction term reads 

(1) ^ " 

4/2 


2 .NN = TVi / dxivn- (r X 7 r) 7 °iV , 


Hf^NN = / dxiV 


^(9*77 ■ (r X 77 ) 7 * + Cl ^^ 77 ^ - -^<9*77 ■ (9*77 -h 

4// fi 


n • n -h ■ ((9*77 X (9j77) 

2/2 


N 


hSnn = / dxiV 


(B9) 


(BIO) 


di + d 2 + ■ 2 ^ 0 A I'-n \ 0 

2- 7 ^ -(r X n) ■ ((9*(9*77 -h 777277 ) 7*^-4 —^(11 x 77 ) ■ tY 

J 7T J n 


+2 


+ ‘'7r X n). (a^s.rr + 7 » + n . S.,r N 7 “ 


f2 
J TT 


J TT 


N , 

(Bll) 
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from which the vertex follows as 


r(i) 

' 27rA 
r(2) 


(^'1 |A;ki,ai;k2,a2) = ^ea,a2cTc (wi - U2) , 


(B12) 


/V\ TjC^) 1X1 1 \_ i 2P + fcr X (p'- p) ml 

(A I \ A; ki, Oi; k 2 , 02 ) — - 2 m- ' ^aia2aT'a + 4ci— ()ai,a 2 


■ —ki ■ k 2 5ai,a2 + 


^CT ■ (ki X k2 ) Cq 


2 c 3 

f? 

C 4 

7, 


2 (c 2 + C 3 ) 

f 2 
J TT 


a;ia;2 <5, 


ai5a2 


(B13) 


(A' I 


(3) 

2-kNN 


A;ki,ai;k2,a2) = i{uJi - CO2) 
^ di -I- d2 


di -|- d2 d^ , d^ m? 


ki ■ k2 + 


P 

J TT 

<^15 ~ <^14 


-a;ia;2 + 4 ■ 


P 


(ki X k2) ■ cr(5a,,, 


a 2 


(B 14 ) 


a 7 n 

and vertices in which either or both pions are in the hnal state are obtained from the equation 
above by replacing kj,a;i —> —kj, —Ui. 


3. 3ttNN vertices 


The interaction terms read 


H. 


(1) _ 9a 


3nNN 


dxiV 


a 77^ T ■ diTT H-(4q; — 1)t ■ TT TT ■ diTT 


PpN , 


(B15) 


which leads to the following interaction vertex 


(A^ I I A; ki, ap k 2 , 02 ; k3, 03 ) — 


The corresponding tadpole contribution is 


(A' I hP^ 


3wNN 


igA 

2/^ 

^ai^a2,a3 [(2 

1/2) 

(k2 + k3) + 2 aki] 

“ 1 “^a 2 ^ai ,a 3 

1(2 a 

- 1/2) (k, 

+ k 3 ; 

) + 2 Q;k2] 


“ 1 “ ^az ^ 0.1 ,a 2 

1(2 0 

- 1/2)(k, 

+ k2; 

) + 2 Q:k3] 

■ (B16) 

L is 






• 9A /-I f\ 

= -«^(10a. 

- 1) Joi Ta 

cr ■ k 


(B 17 ) 


where Jqi has been dehned in Eq. (B2). 


4. vertices 

The interaction term reads 

^Snn = J d^N (n„ 77^ + 77^ n„) (r X 77 )^ 7 ° N (B18) 

and the tadpole contribution follows as 

(0 I I ki, ap k 2 , 02 ) = Jqi ^aia 2 C (a;i-a; 2 ). (B19) 
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5. ANN vertices 


The interaction terms read 

HaU = -y/ dxiVr„^i7 .7"iV , (B20) 

^ANN ~ ~ j dxA^ 1^2 mldiQ r ■ Ai7*75 + ■ (9^Fij7*75^ N , (B21) 

from which the vertices follow as 

(A'l<7|A) = fr.[.-^P^.-j^(p'-p)xP7^..pp 

-W^CT-(p'-p) (p'-p) -Aa, (B22) 

(A' I \ \) = 2ml di 6 Ta(T ■ Aa + ^Ta (qq ■ CT - g^cr) ■ A^ , (B23) 


where in Eq. (B22) terms of order have been retained in the expansion of the bilinear 
N^'j^N, since they have been shown to generate signihcant corrections to the single-nucleon 
axial current [33]. 


6. ttNNA vertices 


The interaction terms read 
^Inna— J dxAAo- (r X 7r)7°A , 

tr(2) _ f _ ^ , A _ ^ _ (.A- sy . A. \ . Aj ^ A * 


(B24) 


HInna = / dx a - —(r X tt) ■ x tt) ■ diAj + -^A* ■ diir 

J L JTT 

+■^(<9*77 X r) ■ Aj N , 

Jtt 

H^^nna = /dx a ^ ^2 + de ^ _ a*AV + ^'9jA° ■ (9i7r(T*^7° 

+2 ^ X tt) ■ did^A^j^ 

Jtt Jtt 


(B25) 


+2 *^1 ^ X Ao) ■ ((9*(9i7r + (9*7rVi)7° -h ^^^ <9^77 ■ Aqo-*^ -f ... A , 

/tt Jtt 

(B26) 


where the dots indicate terms which do not contribute in tree-level diagrams of order Q, for 
example 

/ dx A - 2 7*7^ n ■ d,Ak - r ■ {d^A^ x n) 7 ° A , 


^ d, + d2 + d3 f dxAT-(aoAoxn)7°A, 

JTT J 


45 



and i9oAo —)■ —iUgAo is of order Q^, since in onr connting the energy of the external held is 
of order Q^. The interactions in Eqs. (B24)-(B26) lead to the following vertices 


I hI^nna I k, a) = - ^ eabc n , 


I ^Inna I k, a) — 


4/. 

1 


2 m/, 
. C6 


'^abc 'Tb Ac 


P + -^x(p'-p) 

.•C3 


eabc Tb Ac ■ (CT X q) + 2i— k ■ Ac 

4m/^ U 


.C4 




(3) 

■kNNA 


(B27) 


A; k, a) = 


-i-r^abcnAa- (cr X k) , 

JTT 

2 di - dec * 0 „ 1 , , di4 + 2d23_ 40 


A 


-(A^^ X T)aq - k + 




-cr ■ (q X k)A" 


-^(A»xr)„y 

/tt 


(B28) 


(B29) 


7. 2ttNNA vertices 


The interaction term reads 


H 2 JNNA = -J^ dx AAi ■ [(r X tt) X tt] 7*7^ A , 




which leads to the following vertex and tadpole contribntions 


(A I d727rAfA''A I-^5 kl, fll, k 2 , 02 ) < CO (da,ai ^12 “ 1 “ ^a,a2 2 5a, 

,a 2 ) Aa ■ cr , (B31) 


9A 


(B30) 


4/2 


I H^2nNNA I A) = Joi Ta Aa ■ cr . 


9A 


(B32) 


8 . SnNNA vertices 


The interaction term reads 

r(i) 4a 


H. 


SttNNA 


16/3 


^ dx A 77 ^ A^ ■ (r X tt) 7° A , 


(B33) 


from which the tadpole contribntion is obtained as 

r(i) ix.i. 5(4a-l) 


(A'|i^^:WlA;k,a) = 


32/3 


Joi (r X A°)^ 


(B34) 


9. 27r vertices 


The interaction terms read 


^(4) _ 


dx 


mlk , i X m/ (/g + k) 

(^n ■ n + O TT ■ OiTT) H- - -- TT ■ TT 


f2 
J TT 


f2 
J TT 


(B35) 
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from which the vertex is obtained as 


( 0 | |ki,ai;k 2 ,a 2 ) 


5 


a\,a2 


2 ml U 

P 

L J n 


(uJiLJ2 ~ ki ■ k2) + 


2 ml {h + U) 
P 

J TT 


(B36) 


where, as noted earlier, the momentum-conserving ^-function (27r)^(5(ki + k 2 ) and the pion 
held normalization factor \j are understood. Vertices in which one or both pions 
are in the hnal state follow by replacing a;i,kj —)■ —ca*, — k*. Enforcing the 5 function 
requirement ki = — k 2 = k and iX)\= U 2 = uj, the vertex in Eq. (B36) reduces to 


(0 I hH I k, a; -k, a) 


Amlk 2 , 2m^/3 

f2 ^ + f2 

J TV J TV 


(B37) 


Similarly, we hnd 

(k,a|//£>|k,a> = ^^, (B38) 

J TV 

according to the prescription given above. Apart from the factor 1/(2 ta), which is not 
included in the equations above, these vertices are the same as given in Appendix F of 
Ref. [17]. 


10. Itt vertices 


The interaction terms read 


hS = / rfx 


Aa — 1 


(tT ■ n n ■ 77 -|- (9*77 ■ 77 (9*77 ■ 77) 


„ 0 „ \ 8 (a - 1 O A 

+ — (TT^n ■ nTTa -h 77^(9*77 ■ (9 77 )-—^771^77^ 

J TV TV 


(B39) 


which leads to the following vertex 


(o') -L 

(0 I HlJ I ki, ap k 2 , 02 ; ks, 03 ; k4, 04 ) = 

J TV 

X ^ai,a2^a3,a4. [—2Q;(a;i + U 2 + + (ks + k4)^ -|- (wi -|- U2){u:z + 1 ^ 4 )] 

+<5ai,a3<5a2,a4 [— 2 Q;(a;i -|- ^2 + 1^3 -f + (ki + k 3 )^ -|- (wi -|- a;3)(a;2 + ^^ 4 )] 

+<5ai,a4<5a2,a3 [— 2 Q;(a;i -|- a ;2 + ^3 -f 074 )^ -|- ml -|- (ki -|- k4)^ -|- (wi -7 uja){uj2 + Ci;3)] 


(B40) 


and the corresponding tadpole contribution is 

(o') 

(0 I -^47r I kl, Cll, k2, 02 ) = (5ai,a2 ^01 


1 - 10a , , . 20a -3 21 

{uiU2 - kl ■ ka)-—^ m^ , (B41) 


2 // 


4/^ 


and the constant Jqi has been dehned in Eq. (B2). 
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11. -kA vertices 


The interaction terms read 

^S = A/dx(A*-9,7r + A0-n) , (B42) 

y dx A* ■ 9,77 , (B43) 

from which the vertices are obtained as 

{0|ffS|k,a)=»/,(k.A„-wA2) , (B44) 

{0|i42|k,o) = 2i^k.A„. (B45) 

/tt 

12 . ‘i'lrA vertices 

The interaction terms read 

= 7 ^ / dx 2(1-2 a) A* • 77 77 ■ diTT - (2 a + 1) A* ■ 9^77 77 ■ 77 

+2 (a - 1/2) A^TTftllaTrft + 2a A° ( 7 ra 77 ■ n + n ■ 77 7 ra) , (B46) 

which lead to the following vertices 

(0 I I ki, ai, ka, Oa; kg, og) = ^ 6 a 2 ,a 3 ■ [(2 a - 1 ) q - 2 kg] 

+( 5 ai,a 3 Aaa ■ [(2 a - 1 ) q - 2 ks] 

+( 5 ai,a 2 Aag ■ [(2 a - 1 ) q - 2 kg] 

—< 5 ^ 2,03 [2 a (cui + a ;2 + cug) — cui] 

— Sai,a 3 [2 a (cui + a ;2 + tUg) — C<; 2 ] 

~^a 3 ,a 2 ^^a 3 ^ (cui + a ;2 + ^ 3 ) ~ ^s] 5 (B47) 

where in the hrst three lines use has been made of the ^-function (27r)^(5(ki + k 2 + kg + q). 
The tadpole contribution is found to be 

(0| |k,a) =(5a + 1/2) Aa ■ k+(5a - 3/2)A°a; . (B48) 

Appendix C: Contact terms at order Q 

The weak-interaction potential = A° ps „ — A^ ■ jg „ is parity {V) and time-reversal (T) 

invariant, which implies that —t -p^^a and js^a —t j 5 ,a, and p^^a —t (-)“+^ P 5 ,a and 

'j- 

j 5 ,a —t (—)“j 5 ,a- At Order there is no momentum dependence, and consequently there 
are no contact terms which can be constructed for while two such terms occur for of 
which only one is independent (Fierz identities, see below) and is given in Eq. (5.4). At order 


48 



Q the contact terms in ps „ and must be linear in either kj = p' — pj or Ki = (p- + Pi) /2 
with i = 1 and 2 . None can be constructed for A complete, but non minimal, set of 
hermitian operators for the axial charge is the following: 

Ol = (ri^a + r 2 ,a) + <T 2 ) ■ (Ki + K2) , 

O2 = (ri,a + T2,a) (cTi - <72) ■ (Ki - K2) , 

O3 = i in,a + T-2,a) (cTi X <72) ■ (ki - ks) , 

O4 = {ri^a — T2,a) ((^1 ~ (^2) ' (Ki + K2) , 

O5 = (ri^a — T2,a) 1 + 2) ' (Ki — K2) , 

Oe = i (ri,a - r 2 ,a) (cti x <72) ■ (ki + k 2 ) , 

O7 = i (ti X T2)„ (cTi - <72) ■ (ki + k2) , 

Os = f (ti X T2)„ (cTi + 0-2) ■ (ki - k2) , 

Og = (Ti X T2)„ (cTi X 0-2) ■ (Ki + K2) . 

The antisymmetry of initial and hnal two-nucleon states requires 

Oi = -prp^p^P^ceQ, ^ 

where is the space exchange operator, and P”" and P'^ are the spin and isospin exchange 

operators with = (1 -|- < 7 i ■ <72) /2 and similarly for P’". Exchange of the hnal momenta 
of the two nucleons p'^ ^ P2 leads to 

pspace(k^ + k 2 ) = ki + k 2 , P«P'^“(ki - k 2 ) = 2 (K2 - Ki) , (C 2 ) 

P^P"ce(Ki + K2) = Ki + K2 , p«P--(Ki - K2) = (k2 - ki) /2 , (C 3 ) 

while spin exchange implies 

P'^ (cTi + <72) = <71 + <72 , P”" (< 7 i - <72) = f (< 7 i X <72) , P'^ (< 7 i X <72) = -i (< 7 i - <72) , 

(C 4 ) 

and similar relations follow under isospin exchange. The following (Fierz) identities are 
obtained from Eq. (Cl): 

O2 = O3/2 , O4 = O9 , O5 = Os/2 , Oq = —O'l , (C 5 ) 

while Ol is required to vanish. Hence only 4 of the above 9 operators are independent, and 
a convenient set is 

Ol = ^Or — Os^ /2 , 02 = ^07 -|- Os^ /2 , O3 = (^Oq — Os^ /2 , 04 = 04. (C6) 

We note that Oi and O3 have the same operator structures associated with the divergent 
parts of the loop diagrams. 


Appendix D: Regularized loop contributions to 


MPE 

a 


The regularized contributions of diagrams in Fig. 4 read: 




dz 


<7iM(A;2, z) k2 <7i ■ k2 




M{k2,z) 


(Dl) 
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jS(e4) = 


9a 


jv(®5) = 


kfzz 

Q4^f4 '^'^[M{ki,z) 

rl 


9a 


r 2 ,a cr 2 I dz 

q 


+ 3 M{ki, z) 


dz 


r2,a 0-2 ■ (ki - ks 


1287r/4 g 2 + ^2 
- (ti X Ts)^ (cTi X 0 - 2 ) ■ ki M{ki, z) 


kfzz 

M{ki,z) 


+ 3 M{ki, z) 


jS(e 8 ) = 


9a 


64:71 f* 

/c| 

CTi 

2 



■ 

/ dz 

^ 2 ,a 

Jo 

_ 


5 (TiM{k 2 , z) + —(Ti-k 2 


klizzf ^ l-lzz 


X 


, lzz{l-8zz) 
M{k 2 ,z) ^4 M{k 2 ,zf 
1 ki (2z — 1)^ 


M{k 2 ,zf M{k 2 ,z) 

+ 2 Tin (<^2 X k2) X k2 


jS(elO) = 


[4M(A;2,x) 48 M{k2,zf 

rl 


9a 


dz 


(2 T2,a - Ti^a) 


jS(el5) = 


128 7r/4g2 + 

+ (ti X T2)„ M(/c 2 , z) {(Ti X (T 2 ) ■ k2 

kf zz 


ki 


M{k 2 ,z) 


+ 3Ad{k 2 , z') 


CTi ■ k 2 


9a 


128 7r/4 



■ 

/ dz 

Z~2,a 

Jo 

_ 


M{ki,z) 


+ 3 M{ki,z) 


(k2 — 3 ki 


+4 (ti X T 2 )„ (cTi X ki) M{ki, z) 


0-2 ■ k 2 


u. 


2 ’ 


jS(el 6 ) = 


9a 


q 



■ 

/ dz 

'T2,a 

Jo 

_ 


128 7r/4g2 + ^2 

+3k1 + 3q^ — 20 k'lz'z) + 


— IQ M{ki,zY + M{ki, z){lb ml + 11/c' 
kl zz 


M{ki, z 

-2 (ti X T2)„ (cTi X kl) ■ (k2 + q) M{ki, z) 


(5 ml + kl + q'^ + 3 kl — 2 klzz) 

<T 2 ■ k 2 


2 ’ 

u 1 


‘ ) 22 ^2 

■(1)on'i ^-4 ^71 ^ 

jU(e20) = -^^^C'Tri,aO-2 , 


32 TT g2 + i„2 ^2 


87r/2 

where M{k,z) and ^ have been dehned in Eqs. (5.32) and (5.33). The contributions 
sponding to diagrams e2, e9, and e21 easily follow from those for el, e8, and e20. 
The loop functions Wi and Zi introduced in Eqs. (7.4) and (7.5) are dehned as 


Wi{k) = / dz 


W 2 {k) = / dz 

Jo 


(1 - 5^y M{k,z) - 


9 a 


10;^^- 1 ^ zz{l -8zz) 
M(k,z) ^ 4M(k,z)3 


9 AZzk‘^ ^zz{7g1i + 2 )- g\ 


M{k,z) 


2M{k, z) 


(D2) 

(D3) 

(D4) 

(D5) 

(D6) 

2 

L 

(D7) 

(D8) 

(D9) 

corre- 

(DIO) 

(Dll) 
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w^{k) 

Zi{k) 

^2(k) 


dz 


k^iz-zY ^ 1 


2 Jo [12 M{k,zY M{k,z) 


= dz 


zzk'^ 


0 [M{k,z) 


+ 3 M{k, z) 


(D12) 

(D13) 


= dz 


-20zzkY + 


4ml — 10 M{k, z)^ + M{k, z) (15 + 14 — 6 q • k + 6 

zzk‘^ 


M{k,z] 


(5 + 4A;^ + 2g^ — 2q-k — 2A;^;2 2 ;) 


Zsik) = / dzM{k,z) . 


(D14) 

(D15) 


Appendix E: Counter-terms to order 

Having made the replacements in Eqs. (6.1)-(6.5), the bare Lagrangian C can be rewritten 
in terms of the renormalized helds and physical masses as 


C = C + 6CY 


(El) 


where is the same as in Eq. (2.4) but now in terms of renormalized helds and masses, 
and SC^ includes the set of counter-terms 

5C=5mWN^ + 5ZnN" - nf) + dZ^W [r°’"(0)9o< + + ^"(1)] iV" 


+5Z^N 


dm 


+ 


[r°’^( 0 ) + 5r^(0)] + [Ar(0)/2 + 5Ar(0)] d^TTl + dAYl) 

dZ, ^ 




2 
-^r 2 


Gb + 5Gb d%l + I GA + 5Gb 




-dZ^UA>:{F:j2 + dF:,)d,7;l, 


(E2) 


where rb’'(0), A^’^(O) and A''(l) are the held combinations dehned in Eqs. (A54), (A55) 
and (A56) expressed in terms of renormalized helds and physical masses. The remaining 
quantities are given by 


5r°^(o) = 


So; — 1 




-77 -77 (^T X 77 j 7 , 


5A*/(0) = [2 a 77’’ ■ 77’’ Ta -f (4 a - 1) T ■ 77 ^ 77 ) 


. 7*75 , 


3 a 


5A^(1) = _ )^l - _ 77’- ■ 77’-J (r X 77^) ■ Ao7 

+ 1^ [ ^ ^ b*75 , 

^ 2 O' 1 — 4-0 

jT/or ^ r i 

f2 ^a^b 5 

J n J TV 




(E3) 

(E4) 

(E5) 

(E6) 

(E7) 
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2 a + 1 




It is convenient to define 

5u = S:,+dz,{a:,+dG:,) , 

GL = gU + 2 rTN's^t , 

f’.'i, = + '5^. (h/2 + SF'J , 

= h:, +( j /:^+ 5h:^ , 
rj' = + ^z„rj'( 0 ) + sz, [r^''(o) + ffj'(o)] , 

A-' = Af + iZ„Ai’'(0) + SZ^ [Ai'(0)/2 + Mf(0)] , 
A' = A’’ + «„A’'(1) + ^Z,M’'(1) , 

which then leads to the Lagrangian as given in Eq. (6.6). 
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(E9) 
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Erratum: Nuclear Axial Currents in Chiral Effective Field 
Theory [Phys. Rev. C 93, 015501 (2016)] 


A. Baroni, L. Girlanda, S. Pastore, R. Schiavilla, and M. Viviani 


PACS numbers(s): 21.45.-v, 23.40-s 


Equations (7.2), (D4), (DIO), (Dll), and (D12) of the original article [A. Baroni, L. 
Girlanda, S. Pastore, R. Schiavilla, and M. Viviani, Phys. Rev. 93, 015501 (2016)] are 
erroneous. The correct equations read, respectively: 


• OPE _ 9a 

J5,a 


+ { Ci + 


2C3 + 1 T-2,aka + (Ti X T 2 ) 


128 TT fl 
1 I g\m. 


if. C6 -f 1 

Ki- 7 — CTi X q 


2 m 


Am 


4 m 128 TT 


CTi X k2 


<72 ■ k2 ^ , 




(7.2) 


= 


9a 


647r/4 

/Cn 

<71 

2 



' 

/ dz 

^2,a 

Jo 

_ 

9zz- 

1 


5<7iM(A;2,x) -h -^<7 i ■ k2 


kl{zzY ^ l-7zz 


X 


klizzf 

M{k 2 ,z) M{k 2 ,zY 
1 kl {2z — 1)^ 


M{k2,zY M{k2,z 

2ri„ (cr 2 X ks) x k 2 


+ 


[AM{k2,z) 48 M{k2,z) 


(D4) 


Wi{k) = dz 
Jo 


{1-5 9 \) M{k,z) - ^ 


9zz-l k^izz'^'^ 


\_M{k,z) M{k,zY 


(DIO) 


W 2 ik) = dz 

Jo 


g\{zzYk^ ^zz{7gl + 2)-g\ 


2M{k,zY 2M{k,z) 


(Dll) 


Ws{k) 


- k dz k^z-zf 1 

2 Jo [12 M{k, zY M{k,z) 


(D12) 


In particular, referring to the original article, Eq. (7.2) had the signs of the two terms 
proportional to g\ wrong , and Eqs. (D4), (DIO), (Dll), and (D12) had algebraic errors 
in going from the formal expressions of Eqs. (5.11), and (7.5) of the original article to the 
regularized expressions. These corrections do not affect any of the conclusions of the article. 
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